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ABSTRACT 

Starting in late 2013, the eROSITA telescope will survey the X-ray sky with unprecedented 
sensitivity. Assuming a detection limit of 50 photons in the (0.5-2.0) keV energy band with a 
typical exposure time of 1.6 ks, we predict that eROSITA will detect ~ 9.3 x 10"' clusters 
of galaxies more massive than 5 x 10^^h~^ Mq, with the currently planned all-sky survey. 
Their median redshift will be z ~ 0.35. We perform a Fisher-matrix analysis to forecast the 
constraining power of eROSITA on the ACDM cosmology and, simultaneously, on the X-ray 
scaling relations for galaxy clusters. Special attention is devoted to the possibility of detect- 
ing primordial non-Gaussianity. We consider two experimental probes: the number counts 
and the angular clustering of a photon-count limited sample of clusters. We discuss how the 
cluster sample should be split to optimize the analysis and we show that redshift information 
of the individual clusters is vital to break the strong degeneracies among the model param- 
eters. For example, performing a "tomographic" analysis based on photometric -redshift esti- 
mates and combining 1- and 2-point statistics will give marginal l-a errors of Aas — 0.036 
and AOni — 0.012 without priors, and improve the current estimates on the slope of the 
luminosity-mass relation by a factor of 3. Regarding primordial non-Gaussianity, eROSITA 
clusters alone will give A/nl — 9, 36, 144 for the local, orthogonal and equilateral model, 
respectively. Measuring redshifts with spectroscopic accuracy would further tighten the con- 
straints by nearly 40 per cent (barring /nl which displays smaller improvements). Finally, 
combining eROSITA data with the analysis of temperature anisotropics in the cosmic mi- 
crowave background by the Planck satellite should give sensational constraints on both the 
cosmology and the properties of the intracluster medium. 

Key words: cosmology: cosmological parameters, large-scale structure, early Universe - 
galaxies: clusters - X-rays: galaxies: clusters. 



1 INTRODUCTION 

eROSITA^ (Predehl 2006; Predehl et al. 2010) is the primary sci- 
ence instrument onboard the Spectrum Roentgen-Gamma (SRG) 
satellite^, a fully funded mission with a currently planned launch 
in late 2013. eROSITA will perform an X-ray all-sky survey with 
a sensitivity ~30 times better than ROSAT. Assuming a detection 
limit of 50 photons in the (0.5-2.0) keV energy band and with an 
average exposure time of ~ 1.6 ks, we expect that eROSITA will 
be able to detect all clusters of galaxies in the observable Universe 
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with masses higher than ~3xl0^''/i ^ AIq . This makes it an ideal 
probe of cosmology. 

The evolution in the number density of massive galaxy clus- 
ters as well as their clustering properties strongly depend on the 
cosmological parameters. Cluster number counts trace the nor- 
malization and the growth of linear density perturbations (e.g 
Allen, Evrard & Mantz 2011). The mass function of clusters de- 
tected in X-rays by ROSAT and re-observed with Chandra has 
been recently used to delve into the mystery of cosmic accelera- 
tion and thus constrain the equation-of-state parameter of dark en- 
ergy (Mantz et al. 2010b; Vikhlinin et al. 2009b). These studies are 
based on samples containing a few tens of galaxy clusters. 

The advent of eROSITA will produce cluster catalogs with 
~ 10^ objects. Photometric redshifts of the X-ray-detected clus- 
ters will become available thanks to a series of complementary 
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multi-band optical surveys (e.g. PanSTARRS^, DES", and LSST''), 
while planned massive spectroscopic surveys are designed with 
eROSITA follow-up as science driver, e.g. 4M0ST''. All this will 
vastly improve upon current cosmological constraints. eROSITA 
will likely be the first Stage IV dark-energy experiment (as de- 
scribed in the US Dark Energy Task Force report, Albrecht et al. 
2006) to be realized. 

At the same time, better data will allow us to relax some 
model assumptions. An interesting example is given by the statis- 
tical properties of the primordial density perturbations that seeded 
structure formation in the Universe. They are often modelled with 
a Gaussian random field. In fact, the simplest inflationary models 
predict negligible - i.e. well below the current detection limit - lev- 
els of primordial non-Gaussianity. Nevertheless other theories for 
the generation of the primordial density fluctuations (both infla- 
tionary and not) can produce measurable levels of primordial non- 
Gaussianity: e.g. multiple-scalar-field scenarios, curvaton models, 
ghost inflation, topological defects (see Bartolo et al. 2004; Chen 
2010, and references therein for a complete review). Quantifying 
the primordial non-Gaussian signal has thus emerged as one of the 
most powerful tools to test and discriminate among competing sce- 
narios for the generation of primordial density fluctuations. 

Primordial non-Gaussianity is expected to leave imprints on 
the large-scale structure of the evolved universe, ranging from the 
abundance of collapsed massive objects to their spatial clustering 
properties (see Desjacques & Seljak 2010, and references therein 
for an updated review). Clusters of galaxies as tracers of the un- 
derlying dark-matter (DM) haloes are thus an optimal target for the 
search of primordial non-Gaussianity. 

Regrettably, theoretical models provide robust predictions for 
the halo abundance and clustering only in terms of DM masses. 
Therefore some knowledge of the relation between cluster observ- 
ables and the corresponding halo masses is required to constrain 
cosmological parameters. In X-ray astronomy, cluster masses can 
be determined exploiting a series of scaling relations (luminosity- 
mass, temperature-mass, gas mass-total mass), whose shape and 
magnitude can be inferred by dedicated observations or by invoking 
assumptions about the properties of the intracluster medium (ICM). 
Once the functional form of these relations is known, large cluster 
samples can be used to derive self-consistent constraints on cosmo- 
logical and scaling-relation parameters from the same data set (the 
so-called self-calibration approach, Majumdar & Mohr 2003; Hu 
2003). This should also provide some insight into the gas physics 
of the intracluster medium. Note, however, that dedicated observa- 
tions will be anyway necessary to test the range of validity of the 
assumed scaling relations and provide an independent test of the 
retrieved parameters. 

The goal of this paper is to assess the power of 
eROSITA in the simultaneous determination of cosmological 
and scaling-relation parameters, also accounting for primordial 
non-Gaussianity. Forecasts are made using the Fisher informa- 
tion matrix for the measurements of the abundance and two- 
point spatial clustering of X-ray selected clusters. A series of 
papers have recently addressed similar topics. Oguri (2009) and 
Cunha, Huterer & Dore (2010) discussed optical surveys as the 
HSC'', DES, and LSST Similarly, Sartoris et al. (2010) focussed 
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on the envisaged X-ray satellite WFXT*. More qualitative analyses 
have been carried out by Fedeli, Moscardini & Matarrese (2009), 
Roncarelli et al. (2010), and Fedeli et al. (201 1) for the eROSITA, 
the South-Pole-Telescope, and the Euclid surveys, yet not provid- 
ing detailed forecasts, or with no assessment of the effect of the 
uncertainties in the scaling relations, or only considering a small 
subset of cosmological parameters. Our work differentiates from 
the previous ones in the sense that it is meticulously tailored around 
eROSITA and the actual way observations will be taken. Galaxy- 
cluster number counts and power spectra are calculated in terms of 
the raw photon counts that will be detected at the telescope. This 
is done by adopting observationally motivated X-ray scaling re- 
lations, taking into account the spectral energy distribution of the 
ICM emission, the photoelectric absorption suffered by the photons 
along the line-of-sight, and the expected telescope response. 

The plan of the paper is as follows. In Section 2, we ex- 
plain how we compute the expected abundances and projected an- 
gular clustering of DM haloes in the presence of primordial non- 
Gaussianity and in Section 3 we define our fiducial model. In Sec- 
tion 4, we focus on the transformation between actual halo masses 
and observed photon counts. The eROSITA cluster survey is de- 
scribed in detail in Section 5, while the experiments we consider are 
discussed in Section 6. In Section 7, we describe how we compute 
the Fisher information matrix while, in Section 8, we derive the 
constraints on the cosmological and scaling-relation parameters for 
an all-sky survey with eROSITA, giving special attention to the 
primordial non-Gaussianity parameter of the local type, /nl^'. We 
consider two experiments: the measurement of cluster abundances 
and of the corresponding angular power spectrum. For both probes, 
we distinguish the cases in which redshift information on the indi- 
vidual clusters is available or not. Eventually we combine the dif- 
ferent eROSITA measurements and also discuss how the model 
constraints are consolidated by using the (currently forecasted) re- 
sults of the Planck satellite as a prior probability distribution. Table 
3 summarizes our main findings. We investigate the effects of mod- 
ifying the survey strategies and changing the observational cuts in 
Section 9, while in Section 10 we extend our analysis to other inter- 
esting models of primordial non-Gaussianity, complementary to the 
local type: the so-called orthogonal and equilateral shapes. A thor- 
ough discussion of the assumptions in our analysis is presented in 
Section 1 1 where we quantify the effects of changing the fiducial 
model and adding extra parameters. Our conclusions are summa- 
rized in Section 12. 



2 MASS FUNCTION AND BIAS OF DARK-MATTER 
HALOES 

2.1 Halo definition 

The densest regions of the cosmic matter distribution are com- 
monly dubbed "dark-matter haloes". Their outer boundaries (and 
thus their masses) are somewhat arbitrary depending on the 
adopted definition. In numerical simulations haloes are most com- 
monly identified using either the Friends-Of-Friends (EOF) or 
the Spherical-Overdensity (SO) algorithms (Einasto et al. 1984; 
Davis et al. 1985). While the EOF method approximately tracks 
isodensity contours (Efstathiou et al. 1985; More et al. 2011), the 
SO one fixes the mean density of the structure within a sphere 
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grown around the halo centre. We define Ma to be the mass en- 
closed in a sphere whose mean density is A-times the critical den- 
sity of the universe at the time the halo is considered: 



Ma = g7r(Apc 



(1) 



This definition is widespread among observers and slightly dif- 
fers from what is commonly adopted by the numerical community 
where the background density pm is used instead of the critical 
one. Thus, for a given population of SO haloes, the values that 
the parameter A assumes in the two approaches differ by a factor 

From an observational point of view, it is standard to define 
cluster masses adopting the SO algorithm. In fact, for X-ray obser- 
vations of galaxy clusters, better measurements of the gas mass and 
temperature can be achieved in regions with high-density contrast 
where structures are much brighter and relatively relaxed with re- 
spect to the outer regions (and also the assumption of hydrostatic 
equilibrium is more accurate). On the other hand, observational re- 
sults should be compared with numerical simulations where the fi- 
nite resolution and force softening might create problems for very 
high values of A. For these reasons, the best compromise is to de- 
fine the DM haloes in terms of M500 . 

2.2 Halo mass function 

The halo mass function, dn/dM[M, z), gives the halo abundance 
per unit volume and per unit mass as a function of redshift. It is 
conveniently described in terms of the function / defined as 



— (M, z) = .f{-,z)^ ^ 



(2) 



where pm,o denotes the mean background matter density today, 
and a^{M, z) is the variance of the smoothed linear density field, 

{SUz)): 

g\M, z) = — Pii^ik, z) W\k, M) dfc, (3) 

271-^ Jo 

with Piin(fc, z) the corresponding power spectrum and W{k, M) 
a window function with mass resolution M; here we use a top-hat 
filter in real space, which in Fourier space reads 

sin(a;) — a;cos(a;) 



W{k,M) = 3 ■ 



(4) 



where x = kR and M = ^-KpmR'^. The function / is calibrated 
against high-resolution numerical simulations and different param- 
eterization are available in the literature. Following Tinker et al. 
(2008), we assume that the function / weakly depends on z and 
Ab, where At defines the mean overdensity of SO haloes as in 
Eq. (1) but wrt the evolving mean background density of the uni- 
verse. For Gaussian initial conditions, we thus write 



fT{a,z,Ab)^A 



1 



exp 



(5) 



with 

A{At,z) 
a(Ab, z) 
b{Ab,z) 
c{Ab,z) 



Ar\Ab){l + z) 



TMF 

2 

Co 



ao^^^-(A,)(l + ^)'"^ 
^r"(A,)(l + z) 

TMF/ 



,TMF 



0.75 



logio(A,/75) 



(6) 



where A™^ = 0.14 and a™^ = 0.06. The zero subscripts in- 
dicate the values of the parameters obtained at z = and listed in 
Table 2 of Tinker et al. (2008). Note that because of the different 
definition of SO-halo masses, an interpolation of the best-fit pa- 
rameters of Table 2 in Tinker et al. (2008) is required to compute 
these coefficients. For Gaussian initial conditions, this formula for 
the DM halo mass function reproduces to 5 per cent level of accu- 
racy the abundance of SO haloes in N-body simulations for differ- 
ent choices of the standard cosmological parameters and for a wide 
range of masses and redshifts which encompasses what is probed 
by eROSITA. 



2.3 Non-Gaussian corrections 

We consider here extensions to the standard cosmological model 
where the primordial curvature perturbations (more precisely the 
Bardeen's potential) that generate the large-scale structure are not 
Gaussian. In the synchronous gauge and for sub-horizon pertur- 
bations in the mass density (or in the gravitational potential), 
these models are characterized by a non-vanishing bispectrum, 
B{ki,k2,k^) 7^ 0. The functional form of the bispectrum with 
respect to the triplet of wavenumbers K = {k\, k2, kg) describes 
the "shape" of primordial non-Gaussianity while its overall ampli- 
tude (at fixed K) defines the "strength" of the non-Gaussian signal. 
The latter is generally quantified in terms of the non-linearity pa- 
rameter /JjL where the index i refers to a particular non-Gaussian 
shape. As a reference case, we first focus on the "local" model for 
primordial non-Gaussianity, where the linear gravitational potential 
can be written as 



(x) = 5(x) + /Nr'[5^(x) 



(7) 



where g is a Gaussian random field with zero mean and variance . 
In this case, the bispectrum of linear density perturbations reads: 
B{ki,k2,k3) ~ 2/i°L""' [Piin(fci)Pn„(fc2)+ 2 perms.]. Note that 
we normalize /nl^' by imposing the relation (7) at very high red- 
shift which sometimes is referred to as the CMB normalization. We 
postpone to Section 10 the extension of our analysis to other non- 
Gaussian models, such as the equilateral and orthogonal shapes. 

In general, for non-Gaussian initial conditions, we compute 
the halo mass function through a multiplicative correction to the 
Gaussian one: 



/ dn\ 



/ dn\ 



Rng{M, z, /nl), 



(8) 



(here the subscript T, Ai^oo refers to the fact we use Eq. (5) 
and A — 500). Among the possible choices for the non- 
Gaussian correction R]sig{M, z, f^£^^), we adopt the prescrip- 
tion by Lo Verde et al. (2008), which is in rather good agreement 
with numerical simulations (Pillepich, Porciani & Hahn 2010; 
Giannantonio & Porciani 2010): 



R^g{M,zJI,T') 



dSs 

diner 



4-1 



(9) 



Here 5c is the critical density contrast for halo collapse fixed to 
1.686 for SO haloes (Desjacques & Seljak 2010), while S3 denotes 
the standardized third central moment (skewness) of the smoothed 
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density field, namely Ss = {S'm) /a*' 



a 4 



/local 
NL 

(27r2)2 



dk^kl W{kl,M) Plin{kl 



dkikl W(k2,M) Pnn{k2 



d^j.W{ki2,M) 



1 + 2 



Plin(fcn 
Plin(fc2 



(10) 



where ki2 — ki + k2 + 2^k\k2. 

Note that we do not use the direct fitting formu- 
lae for the /nl''' -dependent mass function presented in 
Pillepich, Porciani & Hahn (2010) as they have been obtained for 
FOF haloes and there is no trivial mapping between FOF and SO 
haloes. 



2.4 Halo bias 

We are not only interested in extracting cosmological informa- 
tion from measurements of the abundance of the X-ray clus- 
ters, but also from the angular clustering of the observed objects. 
From numerical and theoretical work, we know that the cluster- 
ing of DM haloes in which the clusters reside is biased relative to 
that of the underlying mass distribution by an amount which de- 
pends on halo mass, redshift, and the scale at which the cluster- 
ing is considered (see e.g. Mo & White 1996; Catelan et al. 1998; 
Smith, Scoccimarro & Sheth 2007). In general, the bias of DM 
haloes can be expressed as 



h{k, M,z) = 



Pi,Uk,M, z) 



P(k,z) 



(11) 



where Phh is the power spectrum of the halo density field and 
P(k, z) is the matter power spectrum. Within Gaussian scenarios, 
the halo bias on large scales (k < 0.05 h Mpc~^) depends on the 
halo mass and redshift but is independent of k (e.g. Tinker et al. 
2010). We compute this linear bias using the peak-background-split 
model (Cole & Kaiser 1989; Mo & White 1996): 



foe = 1 



1 din/ 
5c diner 



(12) 



On the other hand, a strong fc-dependence 
of the bias appears at smaller scales (fc > 
0.1 /iMpc"^), e.g. Smith, Scoccimarro & Sheth (2007), 
Manera, Sheth & Scoccimarro (2010). 

Dalai et al. (2008) have shown that primordial non- 
Gaussianity of the local type introduces a scale-dependent 
correction to the large-scale bias of massive DM haloes. This is 
because large and small-scale density fluctuations are not indepen- 
dent when /nl^' 7^ 0. We thus extend the expression for the halo 
bias to non-Gaussian scenarios following Giannantonio & Porciani 
(2010): 



M, z) = 6g + /nl (6g - 1) 



a(fc, 2) 



(13) 



which is in good agreement with N-body simulations as shown by 
Giannantonio & Porciani (2010) and Desjacques, Seljak & Iliev 
(2009). In Eq. (13), the different factors read T = 35c llll<? 
and a(k, z) = k'^ T(k) D{z) g{0) /g{oo), where T{k) is the linear 
matter transfer function, D{z) is the linear growth factor of the 
density perturbations (normalized to unity today), and g(z) is the 
growth factor of the potential perturbations. Throughout this paper, 
we use the linear matter transfer function T{k) computed using 



the LINGER code by Bertschinger (2001). 



3 FIDUCIAL MODEL 

We consider a flat ACDM universe characterized by 5 standard pa- 
rameters (3 for the background - Q^n, ^b, and h - and 2 for the 
scalar perturbations - ris and us). On top of this we introduce the 
extra parameter /nl^'. We list all the model parameters and the 
choice of their fiducial values in Table 1, where we distinguish 
among cosmological, cluster-physics (see Section 4), and survey 
(see Section 5) parameters. We assume Gaussian perturbations as 
our fiducial case, i.e. /nl^' = 0- We choose the other fiducial 
values by adopting the best-fitting parameters from the combina- 
tion of 5-yr data from the Wilkinson Microwave Anisotropy Probe 
(WMAP), baryonic acoustic oscillations (BAO), and supemovae la 
(SN), see Komatsu et al. (2009). 



4 FROM HALO MASS TO X-RAY PHOTON COUNTS 

eROSITA will perform an X-ray all-sky survey in about four years 
starting in late 2013, from a L2 orbit. The instrument will consist 
of seven identical Wolter-1 mirror modules of 358 mm of diameter, 
each containing 54 nested mirror shells, and of a fast frame-store 
pn-CCD detector. This will result in a 1 deg field of view, an ef- 
fective area of 1500 cm^ at 1.5 keV, and an averaged point spread 
function (PSF) of 25"-30" half energy width (HEW; on-axis: 15" 
HEW). It will operate in the broad energy range (0.5-10) keV. 

In order to convert non-observable quantities such as M500 
into the raw photons counts effectively collected in the detector (77), 
a series of transformations have to be considered^ : 

— [Lx,Tx] — > photon counts = 77 

where Lx denotes the X-ray luminosity and Tx the average tem- 
perature of the intracluster gas. Consequently, the mass function in 
Eq. (2) has to be converted in a raw-count function: 

dn dn 
dMsoo drj ' 

In full generality, given two variables X and Y, the conversion in 
the differential number counts can be performed as follows: 

^(n^/^wp(y*mdx, (14) 

where P{Y\X) denotes the conditional probability of getting Y for 
given X. When the scatter in the Y — X relation is negligible (i.e. 
the probability P is a Dirac-delta distribution), this reduces to 

dn dn dX 
dY " dXdY' 

On the other hand, any function of the variable X can be written in 
terms of Y as follows: 



(15) 



bYiY") = 



Jb{X)j^{X) P{Y'\X)dX 
j ^{X) P{Y-*\X) dX ' 



which, for negligible scatter and when Y — f{X), gives 

bY{Y'') = b{r\Y')). 



(16) 



(17) 



^ Since the conversion between the energy flux and the photon counts de- 
pends on the energy spectrum, to mimic more closely the experimental re- 
ality we prefer not to adopt a flux limit in our analysis. 
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Table 1. Model and survey parameters plus their fiducial values. The parameters which are allowed to vary in the Fisher-matrix analysis are typed in boldface. 
Note that we assume a flat cosmology where Ha = 1 — f^m throughout the statistical analysis. 



Cosmological Parameter 


Description 


Fiducial Value 


Current error 


Reference 




Non-linearity Parameter (Local) 





—9 < /nl^' +111 


Komatsu et al. (2009) 


CT8 


Normalization of P(fc) 


0.817 


±0.026 


Komatsu et al. (2009) 




Dark Matter Fraction 


0.279 


±0.0058 


Komatsu et al. (2009) 


Us 


Spectral index 


0.96 


±0.013 


Komatsu et al. (2009) 


h 


Hubble Constant 


0.701 


±0.013 


Komatsu et al. (2009) 




Baryon Fraction 


0.0462 


±0.0015 


Komatsu et al. (2009) 




Dark Energy Fraction 


0.721 


±0.015 


Komatsu et al. (2009) 


w 


Equation-of-State Parameter (constant) 


-1 


— U.i4 < i + U? < U.iz 


Komatsu et al. (zUUy; 


X-ray Cluster Parameter 


Description 


Fiducial Value 


Current Error 


Reference 


"LM 


LM relation: Slope 


1.61 


±0.14 


Viklilinin et al. (2009a) 


7LM 


LM relation: 2 -dependent Factor 


1.85 


±0.42 


Vikhlinin et al. (2009a) 


PLM 


LIVl relation: Normalization 




itU.Uoj 


viKniinin et ai. yzxjxjya.) 


CLM 


LM relation: Logarithmic Scatter 


0.396 


±0.039 


ViKnlinin et al. (zOOya) 


Q!TM 


TM relation; Slope 




±0.03 


ViKJilinin et al. (zuuya) 


o 

PTM 


TM relation: Normalization 


Q no \/ 1 nl4 n f J, — 1 
O.U2 X iU IMqil 


r, -1-1 w 1 nl4 
itU.ii X iU 


ViKnlinin et al. (zUUya) 


fTM 


TM relation: Logarithmic Scatter 


0.119 


0.03^* 


Ki-avtsov, Vikhlinin & Nagai (2006) 


Put 


LI correlation coemcient 









ZlCM 


Intracluster metallicity 


O.3Z0 




Anders & Grevesse (1989) 


Nh 


Hydrogen column density along los 


3 X 10^° atom/cm^ 




Kalberla et al. (2005) 


Survey Pai'ameter 


Description 


Fiducial Value 




Reference 




X-ray Energy Band 


0.5-2.0 keV 






'?min 


Minimum raw photon count 


50 








Minimum considered mass (Afsoo) 


5 X IO^Mq/i-i 






/sky 


Sky coverage 


0.658 27, 145 deg^ ^ 




Predehl et al. (2010) 


^cxp 


Exposure Time 


1.6 X 10'' s (all-sky survey) 







" WMAP5-I-BA0-I-SN, for the Cosmology sector (68.3 per cent credibility interval (CI), with the exception of Z^^^' for which the 95.4 CI per cent is indicated) 
^ from hydrodynamical simulations, systematic error encompassing variations for different sub-populations of clusters: relaxed or unrelaxed, at low or high 
redshifts 

all-sky survey excising ±20 deg around the Galactic plane 



In the following, we summarize the specific choices we adopted to 
express the cluster number counts and bias in terms of the photon 
counts. 



4.1 Luminosity-mass relation 

To connect cluster masses to X-ray luminosities, we refer to the 
Lx — M500 relation obtained by Vikhlinin et al. (2009a) through 
the observation of two sets of Chandra galaxy clusters with median 
redshifts of about 0.05 and 0.5, in the (0.5-2.0) keV rest-frame band 
and with Lx integrated within 2 Mpc. The luminosity-mass (LM) 
relation is conveniently written such that the normalization is taken 
at the effective mean mass, to minimize correlations among param- 
eters: 

^jiL = {\nLx) = [/3LM + 1.5(a£M- 0.396^)] (18) 
+ QfLM ln(Af50o/3.9 X lO^Mo) 
+ 7LMlni;(z)-0.391n(V0.72)±crLM, 

where Lx is measured in erg s~^ and E{z) = H{z)/Hq. 
The slope olm reads 1.61 ± 0.14, the normalization 
[/3LM + 1.5(cr^M- 0.396^)] = 101.483 ± 0.085, and the 
redshift-dependence factor 7lm = 1-85 ± 0.42 (see Table 1). The 
symbol ctlm on the right hand side indicates the observed scatter 
in InLx at fixed M. The nature of such a scatter is consistent with 



a log-normal distribution: 



P(\aLx\M) = 



^2^ 



exp 



(19) 



where ctlm = 0.396 ±0.039. This corresponds to a (symmetrized) 
relative error of ~ 40 per cent. 

Note that since Vikhlinin et al. (2009a) determined A/soo from 
the relation between Yx (the product of Tx and gas mass) and 
mass, whose evolution was assumed to be self-similar based on re- 
sults from numerical simulations, the redshift evolution of relation 
(18) has not been measured directly. 



4.2 From X-ray luminosity to the number of detected 
photons 

Given Lx in the (0.5-2.0) keV rest-frame energy band and the 
redshift z of the clusters, we can predict the number of photons 
detected by eROSITA in the (0.5-2.0) keV observer-frame energy 
band. To do this, we require further information: (i) the emitted 
spectrum, which depends on the temperature and metal abundance 
of the intracluster gas; (ii) the photoelectric absorption suffered by 
the X-ray photons along their way; (iii) the eROSITA instrumental 
response and effective area; and (iv) the exposure time, Texp. 

(i) We derive the cluster average X-ray temperature using the 
results by Vikhlinin et al. (2009a), where the relation AIsqq-Tx is 
parameterized with a power law (as expected if the temperature of 
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the intracluster medium scales with the depth of the gravitational 
potential). To measure this relation, Vikhlinin et al. (2009a) define 
the cluster temperature as the average spectral temperature inte- 
grated within a given radial range. Inverting the scaling-relation of 
their Table 3 with no pre-fixed slope, we write the temperature- 
mass relation (TM) as 

= (inTx) = axM ln(M5oo//3TM) + otm lni5(2) 

+ ln(5keV) ± cttm (20) 

where Tx is measured in keV, the slope qtm = 0.65 ± 0.03 
and ^TM = (3.02 ± 0.11) x lO^/i^^Af©. The symbol cttm 
is the scatter in InTx at fixed mass, and reads ctlm — 0.119 
(see Table 1). This value is derived from numerical simulations 
(Kravtsov, Vikhlinin & Nagai 2006) and we assume a lognormal 
distribution also for the conditional probability of getting Tx given 
M500, in analogy with Eq. (19). The fit in Eq. (20) holds for re- 
laxed clusters only: to take into account unrelaxed morphologies, 
the estimated mass from Eq. (20) should be multiplied by a fac- 
tor of 1.17. We do not implement this correction. For the metal 
abundance, we choose a value of 0.3 relative to solar abundance 
(Anders & Grevesse 1989), which is typical for nearby clusters, 
with little or no evolution to higher redshifts (Tozzi et al. 2003). 
We use the apec model (Smith et al. 2001) within XSPEC (Arnaud 
1996, version 12.5. In) to model the emitted spectrum (thermal 
bremsstrahlung, recombination emission, as well as bound-bound 
transitions from collisionally ionized, optically thin plasma), for 
given temperature, metal abundance, luminosity, and redshift. 

(ii) The photoelectric absorption can be modeled using the hy- 
drogen column density along the line-of-sight, A'^h, and assuming 
the associated abundance of helium and heavy elements. For this, 
we use the phabs model within XSPEC (with the default cross sec- 
tions), in order to self-consistently use the same metal abundances 
as above. The Galactic TVh varies across the sky. We use recent 
data from 21 cm radio measurements provided by Peter Kalberla 
(the Leiden-Argentine-Bonn, LAB, survey; Kalberla et al. 2005) 
and determine the A^h distribution, excising ±20 deg around the 
Galactic plane. We find that 3 x 10^" atoms cm^^ is a typical value, 
which we adopt throughout. 

(iii) We model the eROSITA response using the latest pre- 
liminary matrices provided by Frank Haberl^". In particular, we 
use the file erosita_iv_7telfov_ff.rsp, which takes into account the 
CCD quantum efficiency and filter transmission, as well as the mir- 
ror area of all seven telescopes, averaged over the field-of-view. 
The response of eROSITA in terms of count rates is shown in 
Fig. 1 for some characteristic redshifts and for a cluster luminos- 
ity of Lx = W'^ergs'^ 

(iv) The exact orbit and scanning strategy of eROSITA are 
currently under investigation. We consider two different versions of 
the expected four-year all-sky exposure map kindly provided by Jan 
Robrade. We exclude ±20 deg around the Galactic plane. We then 
calculate the average exposure in the remaining area, which results 
in 2.4 ks, for both exposure map versions. This exposure assumes 
100 per cent efficiency. For a more realistic estimate, we need to 
multiply 2.4 ks by a factor smaller than 1. This factor is rather un- 
certain at this time. An optimistic value is 0.8, a pessimistic one 
0.8 X 0.67 — 0.54 (the factor 0.67 is supposed to account for soft 
proton flares, which are expected at eROSITA's L2 orbit). There- 
fore, we expect 1.3 ks < Tcxp < 1.9 ks. We assume 1.6 ks as a 
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Figure 1. Characteristics of the eROSITA telescope: eROSITA re- 
sponse (rate of X-ray photons counted in the detector) in the (0.5-2.0) 
keV energy band (obsei'ver frame), for a source of luminosity Lx = 
lO^^erg in the (0.5-2.0) keV rest-frame band, at different redshifts, 
z = 0.26, 0.31, from top to bottom. 



realistic estimate of the average exposure time, unless specifically 
mentioned otherwise. 

Once Lx , Tx and z are known, we are thus able to compute 
the expected photon counts fja^iLx ,Tx , z) within Xspec for our 
fiducial reference cosmology. Since r; depends on the assumed cos- 
mology through the luminosity distance Dl, we scale the photon 
counts according to = rj&dD^.nd/ D\, every time the assumed 
cosmology is varied. We also account for the statistical uncertainty 
in the photon counts by assuming that they follow Poisson statis- 
tics. In brief, the probability distribution of the counts is 



P{ri\M, z)^ j 



dLx dTx P[v\v{Lx , Tx , z)] P{Lx , Tx | M, z) , 

(21) 

where 77 is a Poisson variate of mean and P generically denotes 
any probability density function. It is interesting to consider the 
possibility that there is a correlation between the X-ray tempera- 
ture and luminosity for clusters of a given mass. In fact this would 
reduce the scatter in the distribution of ri at fixed mass. We write 
the joint conditional probability P{Lx ,Tx\M, z) as a bivariate 
lognormal distribution 



P(X|M) 



27r!S|i/: 



■ exp 



--(x-m)^s-^(x- 



(22) 



where the vector X and its covariance matrix read 



X = 



InLx 
InTx 



and S 



„2 

pLT Ctlm Cttm 



Pltctlmcttm 
Cttm 



(23) 

Here /i refers to the mean values of the scaling relations and plt 
is the linear correlation coefficient between the residuals in lumi- 
nosity and temperature. Theoretically it is expected that Lx and 
Tx may be highly correlated (Stanek et al. 2010) or anticorrelated 
(Kravtsov, Vikhlinin & Nagai 2006); on the other hand, data for a 
set of clusters drawn from the ROSAT all-sky survey with Chandra 
follow-up show no evidence of a correlation (Mantz et al. 2010a). If 
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Figure 2. Observable-mass relation and sample selections for the eROSITA galaxy clusters. Top left: effective photon count vs mass relation (t; — M500) 
derived via Eq. (21). The soHd curves indicate the mean relation at 2 ~ 0.01, 0.10, 0.36, and 1.51, from top to bottom, along with a 68 per cent level of 
scatter encompassed by the dashed lines (corresponding to about ± 41 per cent symmetrized relative error in the range of interest). The black solid line marks 
the detection threshold of 50 photons. Top right: mass threshold (M500) corresponding to different photon-count detection limits as a function of redshift. We 
choose 50 photons with an integration time of 1.6 ks as our reference choice. The thick horizontal line indicates the additional cut at Af = 5 X W^'^h~^ Mq 
which we impose to avoid considering clusters that possibly do not obey power-law scaling relations or exhibit a mass-dependent scatter. The colored curves 
show the effect that the current uncertainties on the parameters of the luminosity-mass relation (see legend) would have on the determination of the detection 
limit. Bottom left: threshold luminosity as a function of redshift, for different photon-count detection limits. The colored lines indicates the position of the 
additional possible mass-cuts in the luminosity-redshift plane (see legend). Bottom right: selection criteria in the r) — z plane. The black solid lines correspond 
to the detection thresholds of 50 and 500 photons; the colored curves connect points corresponding to a given mass limit (colors as in the left panel). 



we assume that Lx and Tx are independent variables and naively 
combine the mass-luminosity and mass-temperature relations from 
Vikhlinin et al. (2009a) we find that the (intrinsic) scatter in lumi- 
nosity at fixed Tx should be ~ 0.49 (in the In-ln plane). In case 
of a positive correlation among the variables, this scatter should be 
smaller. Direct measurements of the scatter from the REXCESS 
survey (Pratt et al. 2009) give 0.276 for core-excised clusters and 



0.666 when the cluster cores are included. On the other hand, other 
studies of the 64 brightest clusters in the sky (HIFLUGCS see 
Mittal et al. 2011) suggest a reduced importance of cool cores as 
scatter contributors and an intrinsic scatter in the Lx-Tx relation 
of 0.455 for the whole sample. Estimating the degree of correlation 
by directly comparing these figures is thus inconclusive and better 
data are needed. We have initially attempted to include plt in the 
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set of parameters that we would like to constrain with eROSITA. 
However, we found out that the weak dependence of the photon 
counts, 7], on Tx (see Fig. 1) makes all the observables rather insen- 
sitive to Put which remains unconstrained by the data. To simplify 
the analysis, we thus assume pur = 0. 

The resulting photon counts as a function of cluster mass and 
redshift are shown in the upper-left panel of Fig. 2. This gives an 
effective rj ~ Afsoo relation which scales as M° with a ~ 1.60 — 
1.65 (depending on redshift). The scatter in our observable-mass 
relation at fixed mass is dominated by ctlm but also depends on 
ffTM and the Poisson noise in the photon counts. This corresponds 
to a (symmetrized) relative uncertainty at fixed mass of about 41 
per cent at 2 ~ 0.1; at 2: ~ 1 this ranges between 66 and 41 per 
cent for objects with M > 5 x lO^/i'^M©. 



5 THE eROSITA EXPERIMENT 

For our purposes, we characterize the eROSITA survey using three 
numbers: the fraction of the sky covered by the experiment /sky, 
the exposure time Tcxp, and the detection threshold in raw photon 
counts rimin (see Table 1). We take as a reference an all-sky survey 
covering more than 27,000 deg^ and with an average exposure time 
of ~ 1.6 ks; moreover, we limit our analysis to the (0.5-2.0) keV 
X-ray energy band, such that all the figures in this work have to 
be considered in that band. This all-sky survey should exquisitely 
constrain model parameters, since it will simultaneously provide 
robust statistics (~ 10^ galaxy clusters detected with more than 
50 photons) and sample very large spatial scales where primor- 
dial non-Gaussianity is expected to leave its strongest imprints 
(k < 0.01^ Mpc~^). As a reasonable source-detection threshold 
we choose r; = 50 counts as the default. This number is motivated 
by the experience from the RO SAT All-Sky Survey (RASS), where 
the count limit has typically been set at (10-30) photons in scien- 
tific analyses of cluster catalogs, e.g., Schuecker et al. (2001). In 
the future, this limit will be tested against detailed simulations of 
the eROSITA survey. While 50 source counts are sufficient for a 
luminosity determination, one could measure the intracluster gas 
mass already with about 500 source counts. The gas mass is a bet- 
ter proxy of the total gravitational mass than the luminosity, so 
when applying the 500 count limit one could in practice assume a 
smaller scatter in the observable-mass relation. While in this Sec- 
tion we discuss both cluster samples determined using these dif- 
ferent thresholds, the whole statistical analysis of Sections 7 and 8 
will be performed with the cut of 50 photons. 

Low-redshift objects which are fainter than lO*'^ erg 
could in principle be detected with more than 50 photon 
counts with an average exposure time of 1.6 ks (Fig. 2, 
bottom-left panel). This corresponds to masses smaller than ~ 
10^^ Mq at z < 0.1 (same Figure, upper-left and upper- 
right panels). Although scaling-relation studies which extend to- 
wards masses of ~ lO"/i"^M0 are available (Sun et al. 2009; 
Eckmiller, Hudson & Reiprich 2011), the adopted prescriptions 
from Vikhlinin et al. (2009a) are obtained from a cluster sam- 
ple limited at about 7 x 10^^ Mq. Moreover, the scatter in 
the observed relations increases for low mass clusters and groups 
(Eckmiller, Hudson & Reiprich 201 1): this would necessitate to in- 
troduce a mass dependent scatter in our analysis. For these reasons, 
in addition to the photon-cut of 50 photons assumed for cluster de- 
tection, we add a redshift-dependent luminosity cut, which is effec- 
tive in the low redshift part (z < 0.25) to remove low-mass systems 
(M < 5 X lO^^ft-^Mo): see Fig. 2, upper-right and bottom-left 



panels. In practice, we pick our fiducial cosmology and our fidu- 
cial luminosity-mass relation as a function of 2 and we look for a 
redshift-dependent luminosity cut that corresponds to the desired 
minimum mass: this, in turn, results into a redshift-dependent 77 cut 
(Fig. 2, bottom-right panel) which is model independent and which 
can be easily applied to real data (while a true-mass cut cannot), if 
redshift information is available Note that neither a luminosity 
cut nor a mass cut are in practice applied here: low-mass objects 
are removed solely referring to the unique curve in the 77 — 2 plane 
which corresponds to the chosen mass threshold in the fiducial 
model (Fig. 2, bottom-right panel): such cut is kept fixed through- 
out the statistical analysis independently of the assumed values for 
the cosmology and scaling-relation parameters (see Section 9.5). 

In Fig. 2, we compare a series of photon-count thresholds 
(■q —50 and 500 photons with Tcxp = 1.6 ks, and 77 =50 photons 
with Tcxp ~ 1.6, 3.0, 7.5 ks) and additional luminosity cuts (corre- 
sponding to masses of~ 1,5, 10 X lO^^ft^M/0), intheM5oo-2 
plane (upper-right panel), in the Lx — 2 plane (bottom-left panel), 
and in the 77 — 2 plane (bottom-right panel). Fig. 2 shows that, away 
from the Galactic plane, eROSITA will detect all clusters more lu- 
minous than ~ 5 X 10** erg up to 2 ~ 1.5, namely all massive 
(M > 3 X 10^* Mq) clusters in the observable Universe. 

In the left panel of Fig. 3, we show the redshift distribution of 
the population of clusters detected by eROSITA (see Section 5 and 
Eq. (24) for details): the solid black curve is the redshift distribution 
of all the galaxy clusters above the cut of 50 photons for Tcxp = 1.6 
ks; the corresponding total number of objects in the whole sky 
(/sky = 0.658) is ~ 1.37 x 10^ objects. In the same panel, the 
dashed curves indicate the redshift distributions of the same pop- 
ulation of clusters once the additional cuts are applied. The inte- 
grated number of objects reduces to ~ 9.32 x 10* and the median 
redshift shifts to ~ 0.35 when objects below 5 x lO^^/i'^M© 
are removed. We summarize our findings in Table 2. To facilitate 
comparisons with other works using different observables than raw 
photon counts, we show in the right panel of Fig. 3 the cumulative 
mass distribution of clusters above our detection threshold of 50 
counts, as function of 2 and M500. 

Finally, in Fig. 4, we compare the selections discussed so far 
with an alternative flux-based selection, often employed in X-ray 
cluster surveys. Two cases are shown: a comparison with a 50- 
photon limit (main box) and with a 500-photon limit (smaller In- 
box). The solid curves show the fraction of objects which pass a 
threshold both in flux and in photon counts with respect to the 
objects which pass a threshold set only in terms of flux. Mini- 
mum fluxes of 3.2,4.2,4.4 x 10"" erg s"^ cm"^ in the (0.5- 
2.0) keV energy band would result, respectively, in 68, 90, 95 
per cent completeness levels for 77niin = 50. For rjmin = 500, 
68, 90, 95 per cent completeness is reached with flux limits of 
3.5, 4.4, 4.6 X 10~^^ erg cm"^, respectively. 



6 OBSERVABLES 

The measurements that we consider are galaxy-cluster abundances 
and angular clustering. In this section we study their sensitivity to 
the actual values assumed by the model parameters. 

Fig. 5 shows how the eROSITA number counts depend on 

In order to apply this additional cut redshift measurements ai'e only re- 
quii'ed for z < 0.25 where, most likely, no new clusters will be discovered 
by eROSITA. Therefore this sample selection should be easily feasible in 
practice. 
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Figure 3. Left: Redshift distribution of eROSITA clusters for different minimum masses (color coded as in tlie bottom panels of Fig. 2). Right: Number of 
clusters passing the detection threshold of 50 counts, as a function of redshift and A/500, with sky coverage fraction /gj^y = 0.658. 











-^clusters 
„ = 50, Toxp = 1.6 ks) 


■^median 


all objects 








1.37 X lO'' 


0.25 


M > 1 X 10" 


{h- 






1.31 X 10^ 


0.27 


M > 5 X 10" 


{h- 


1M0) 




9.32 X lO* 


0.35 


M > 1 X IQi* 


(h- 


1M0) 




5.57 X 10"* 


0.46 



Table 2. Number of galaxy clusters detected by eROSITA and their me- 
dian redshift for the all-sky survey, with sky coverage fraction /^i^y = 
0.658. 



the cosmological parameters and on the parameters characterizing 
the physics of the ICM. In the upper panels of Fig. 5, the "count 
function" ^ at the median redshift of the sample is shown for our 
fiducial model (black solid line) and for alternative models where 
the cosmological (left) and ICM (right) parameters are varied (one 
at the time) by the current uncertainty (listed in Table I). 

The redshift distribution of the clusters above the detection 
threshold of jy^in — 50 counts with Tcxp ~ 1.6 ks relates to the 
count function as follows 




10 



10-"" , ,10-"" 

S . (erg s cm ) 



Figure 4. Completeness of a survey limited in flux and photon counts wrt 
a purely flux-limited selection with S > S^nin f™ a 50-photon limit (main 
panel) and 500-photon limit (inset). 



dz 



dN 47r 

2(,>J?min,Z) — . /sky 



f 

■J Vrr 



H{2 



dn 
dr] 



(24) 



where A is the survey area in deg^ , c is the speed of light, and Da 
is the comoving angular diameter distance (which coincides with 
r for flat universes). The total number of clusters detected above 
a certain detection threshold can be calculated by integrating the 
equation above over redshift and multiplying the result by A. 

The bottom panels of Fig. 5 give the relative change of the 
observables wrt the fiducial model. The largest deviations are 
due to changes in erg, Jim, and and all the four parameters 

of the LM relation. All parameters but tis and h generate larger 
uncertainties at higher redshifts than at lower redshifts (not shown). 
Note that changing as within the current uncertainty modifies 
the cluster counts by 20 per cent for objects detected with more 
than ~ 400 photons at the median redshift, or for all the objects 
above z > 0.6. Also primordial non-Gaussianity (with positive 
/nl''') has a stronger impact at higher redshifts and for higher 
photon counts. The effects of the Hubble constant and of the scalar 
spectral index exceed the per cent level only for the brightest 
clusters, e.g. above a few 10^ photon counts around the median 
redshift, which already suggests that binning the eROSITA data 
in rj should better constrain these parameters. Current uncertainties 
in the parameters regulating the temperature-mass relation (not 
shown but discussed in Section 11.3) affect the cluster counts at 
the per cent or sub per cent level thus indicating that eROSITA 
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Figure 5. Number counts of eROSITA clusters. The black solid line refers to our fiducial model while the other lines show the dependence on the cosmological 
parameters (left) and on the parameters of the LM scaling relation (right), varied one at the time of the relative amounts indicated in the legend (variations 
correspond to the current uncertainties on the parameters as indicated in Table 1). 



counts will hardly be able to put interesting constraints on those 
parameters. Finally, we have checked that fitting errors in the 
parameters of the Tinker mass function (not shown but discussed 
in Section 11.1) also give sub per cent modifications. This confirms 
that our theoretical framework is robust. 

We now consider the spatial clustering of the detected clusters. 
Imagine to split the eROSITA sample into bins (based on the clus- 
ter redshift or photon counts). The angular cross spectrum between 
the pair of bins i, j, can be written as 

c\ roci poo 

Ce(i,j) = -/ dfcfcM drlW^(rl,k) je{kri) (25) 
^ Jo Jo 

X / dr2Wj{r2,k) ji{kr2) PNL{k,r\,r2)- 
Jo 

where the functions Wi (r, k) are defined as 

W,(r,fc) = ^^(r)&.(fc,z(r)), (26) 

and Ni is the total number of objects within the i-th bin, ji denotes 
the spherical Bessel function of the first type, and bi[k, z{r)) is 
the effective bias of the cluster population (obtained by averaging 
Eq.(13) over halo mass or photon counts and weighting by the num- 
ber density) in the i-th bin. Here Pnl denotes the cross-spectrum 
between the matter distribution at two different redshifts, zi and 22, 
corresponding to the comoving distances, r'i(zi) and 7-2(22)- 

In Fig. 6 the dependence of the angular power spectrum on 



the cosmological and ICM parameters is shown for galaxy clusters 
with 77 50 and lying within the broad redshift range [0.1-1.5]. 
Note that the largest modifications in the signal are driven by the 
current uncertainties in cts, Jim, /nl^' and 7lm- An accurate deter- 
mination of the cluster power spectrum with eROSITA has there- 
fore the potential to strongly constrain these model parameters. 

To speed-up calculations it is convenient to rewrite Eq. (26) 
using the Limber approximation: 

(27) 

with the Limber weight functions being defined as 

, , 1 dN, f £ + l/2 \ 

V is the comoving volume. We use Eq. (26) for all multipoles with 
£ < 100 and Eq. (27) for larger values of £. A motivation for this 
choice and a critical discussion about the limitations of the Limber 
approximation in the presence of primordial non-Gaussianity will 
be given in Section 11.4. 

6.0.1 The choice ofimin 

Primordial non-Gaussianity of the local type modifies the shape 
of the angular power spectrum of galaxy clusters by introducing 
a scale-dependent bias on very large scales (see Fig. 6). There- 
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Figure 6. Angular power spectrum of eROSITA clusters and its dependence on the cosmological and ICM parameters. The signal refers to clusters detected 
with »7 ^ 50 photons (Texp =1.6k:s) within the redshift range 0.1 < z < 1.5. 



fore, experimental constraints on /nl^' improve by measuring clus- 
tering at larger and larger scales. However, the correlation func- 
tion of biased tracers of the cosmic density field in scenarios with 
/nl''' 7^ formally diverges for every spatial separation if Fourier 
modes down to fc — )■ are considered (e.g. Wands & Slosar 2009). 
This is because the biased density field has infinite variance (as like 
as the gravitational potential which generates the scale-dependent 
bias (see Giannantonio & Porciani 2010). Any observed correla- 
tion function, anyway, will be finite because density fluctuations on 
scales larger than the survey are not observable. Fluctuations are 
always defined with respect to the mean density as measured from 
the same survey, thus forcing their average to zero. The expectation 
value of the observed correlation therefore departs from the under- 
lying one by a constant shift (sometimes known as the "integral 
constraint") which is finite only if /nl^' = 0. In practice, this ef- 
fect is taken into account by including the window function of the 
survey in the calculation of Ci. This regularizes all the integrals. 
For simplicity, here we approximate the influence of the window 
function by using a high-pass filter k > fcmin ~ 10~^/iMpc~^ and 
only considering a minimum order for the spherical harmonics £niin 
so that the derivatives in the Fisher matrix do not depend on fcmin 
(see also Cunha, Huterer & Dore 2010). In what follows, we will 
use fniin = 5 for the eROSITA all-sky survey, and ^min = 7-10 
for the deeper surveys presented in Section 9.1. These choices are 
motivated as follows. First, multipoles with I < IQ can not be 
precisely measured if /sky substantially departs from unity, since 
only a few modes are available in the whole sky. Second, evalu- 
ating Cf at such low multipoles requires knowledge of the halo 



bias at wavenumbers k <C 0.01 h Mpc~ , never probed by N-body 
simulations. Pushing the analysis to £ ~ 3 would imply trusting 
the extrapolation of the non-Gaussian halo bias of Eq. (13) down to 
k ~ 10^* /iMpc"^ where also general-relativity corrections cer- 
tainly become very important (e.g. Yoo 2010; Baldauf et al. 201 1). 

6.0.2 The choice ofima.^ 

We also have to set a minimum angular scale that will be considered 
in the clustering analysis. This is mainly dictated by the limitations 
of the theoretical models for the angular power spectrum of galaxy 
clusters. Angular multipoles and wavenumbers are related through 
the angular diameter distance in a redshift-dependent fashion. The 
choice of the maximum multipole thus corresponds to selecting the 
largest wavenumber that we want to consider, fcmax- Three issues 
have to be taken into account: cluster exclusion effects, dynamical 
non-linearity of density perturbations, and non-linearity of the DM 
halo bias. We compare the characteristic scales of these three ef- 
fects in Fig. 7. Galaxy clusters are extended objects (at z ~ 0.35, 
a spherical 10^^h~^ Mq cluster has a comoving R^oo of about 0.6 
h'^ Mpc) and their spatial separation cannot be much smaller than 
their characteristic size (~ 2_R5oo)- This imprints a sharp drop in 
the amplitude of Ct at small angular scales which is difficult to 
model accurately. The blue dashed curve in Fig. 7 indicates the 
characteristic size (and thus, roughly, the minimum separation) of 
objects with M ~ IO^^/i'^A/q. Regarding the non-linearity of 
the matter-density field, we do not want to rely on approximate 
prescriptions for the matter power spectrum on mildly non-linear 
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Figure 7. Maximum wavenumber which can be robustly modeled in the 
power-spectrum study, as a function of redshift. The blue dashed line refers 
to exclusion effects and indicates the minimum separation allowed for 
(spherical) objects of ~ lQ^'^h~^ Mq. The black solid line marks the on- 
set of the non-linear regime for density perturbations (see main text for 
details). The grey area approximately traces the minimum wavenumber (as 
a function of redshift) at which a scale-dependence in the bias of dark- 
matter haloes of mass ~ 3 X lO" < A/ < 2 X W^'^h'^MQ can be 
detected in the numerical simulations of Pillepich, Porciani & Hahn 2010. 
Note that, for galaxy clusters, the most severe requirement comes from the 
non-linearity of the bias. To overcome this problem, in our study we only 
consider wavenumbers such that k < fcmax = 



scales which could compromise the accuracy of the fits for the cos- 
mological parameters. The onset of the non-linear regime for the 
matter power spectrum can be approximately evaluated by deter- 
mining fcmax(2) such that 



I Punik,z)dk^0.5 . (29) 



The solid black line in Fig. 7 shows the function fcmax(z). Simi- 
larly, we are concerned with the linearity of cluster bias. Eqs. (12) 
and (13) asymptote to a constant value for k — >■ oo. However, 
N-body simulations show that on sufficiently small scales (de- 
pending on redshift and mass), the bias of DM haloes departs 
from a constant thus becoming non-linear. For haloes with mass 
M > 5 X lO"/i"^M0 at 2 ~ 0, the non-linear regime of the 
bias kicks in for wavenumbers above k ~ 0.1 — 0.3 h Mpc~^. 
This threshold decreases with redshift and mass. The grey band 
in Fig. 7 indicates an approximate determination of the onset 
of non-linearities in the bias for DM haloes in the mass range 
(3 - 20) X 10"ft~^A/o, as extracted from the N-body simula- 
tions of Pillepich, Porciani & Hahn (2010). The numerical work by 
Manera, Sheth & Scoccimarro (2010), shows that for objects with 
M > lO^^/i~^Af0 the transition occurs at even slightly smaller 
wavenumbers. In summary, at all redshifts greater than 0.1, the 
linearity of the bias gives the most stringent constraints, followed 
by the linearity of the density field. Based on these arguments, we 
make a conservative choice and adopt a maximum wavenumber of 
fcmax = 0.1 ^ Mpc~^, which corresponds to £max ~ 96 at the me- 



dian redshift of the all-sky survey. We relax this choice in Section 
11.5. 



7 FISHER-INFORMATION FORMALISM 

The Fisher- matrix formalism (Tegmark, Taylor & Heavens 1997) 
is a tool for forecasting how well a future experiment will constrain 
some model parameters. Given a dataset, x, and a set of parameters 
with fiducial values, 0q, we define the Fisher information matrix 

as 



a/3 



(30) 



where the average is taken over an ensemble of realizations of the 
experiment and C = — InL with L{x; 0) being the likelihood of 
the data given a model. For unbiased estimators, (O) — Oq (where 
6 denotes an estimate of the parameter set based on a single realiza- 
tion of the experiment), the inverse of the Fisher matrix (evaluated 
at = Oo) provides an approximation for the covariance matrix of 
the model parameters: 



y^param 



)a0- 



(31) 



Under some weak regularity conditions, the Cramer-Rao inequality 
assures that the marginal errors of the model parameters follow 



1/2 



(32) 



while the conditional errors (i.e. the errors that would be obtained 
keeping all the other parameters fixed) will be larger or equal than 

1/ V Fact- 

For a (multi-variate) Gaussian likelihood function L{x;0), 
the Fisher matrix can be explicitely written as 

Fal3 = iTr|^S"^(/i,Q fl^,/! +/X,/3 ^J.'^,a ) + S"^E,a S"^S,,3 j 

(33) 

where fi — (x), S — {{x — fi){x — fi)'^), and commas denote 
partial derivatives wrt the model parameters. If A4 is the number of 
model parameters and A/" is the number of data points, the Fisher 
matrix is a x M matrix while S is A/" x A/". The trace operator 
in Eq. (33) acts on the A/"— dimensional space of data points. 

If different experiments are independent, the total Fisher ma- 
trix is the sum of individual Fisher matrices. "Adding priors" can 
be regarded here as adding the Fisher matrix of a third experiment. 



7.1 Figure of merit 

In order to compare individual probes, optimize the experiments, 
and quantify their individual performance in constraining model 
parameters, we introduce the (total) figure of merit (FoM), defined 
as 



FoM = logio [det(F' 



-1/2 



(34) 



This is inversely proportional to the volume of the A^ —dimensional 
error hyper-ellipsoid. The higher the FoM, the more suitable an 
experiment is to constrain our selected parameter set. For example, 
halving uniformly all the parameter errors would correspond to an 
increment in the figure of merit of logig (2^ ) , which for Al — 10 
reads ^ 3. Similarly, if one wants to focus on a selected subset of 
parameters, a (partial) FoM can be defined by considering a Fisher 
matrix of lower dimensionality. 
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7.2 Number counts 

For an unclustered point process, the probability of counting Ni 
objects in the i-th bin is given by the Poisson distribution 



p{Ni\fii) = -jj-fi^' e 



(35) 



where the mean fi = (AT) is the average over an ensemble of re- 
alizations. Fluctuations in the counts are larger for point processes 
that display spatial clustering. In this case the covariance matrix of 
the binned counts can be written as (Lima & Hu 2004) 



(36) 



where S denotes the sample covariance encoding information from 
the two-point correlation function, and the second term is the diag- 
onal Poisson noise Mij = {Ni)Sij. Hu & Kravtsov (2003) have 
shown that the sample covariance is subdominant wrt to Poisson 
errors when a cluster survey encompasses a large fraction of the 
sky, covers a broad redshift interval, and high-threshold masses 
are adopted. We therefore neglect it in our calculations as we will 
treat clustering separately as described in the next section. Note 
that this is different from considering a count-in-cells experiment 
where the covariance of the counts is a direct observable (Oguri 
2009; Cunha, Huterer & Dore 2010). For negligible sample covari- 
ance, the Fisher matrix for the cluster counts can be written as 



f: 



O10 



Tr[St,ounts(Mia /^"^i/s m'^iq )] = 
V), 

bins 

dOa (Ni) 86 13 i fiducial model 



2 

{N),lM-\N) 



./3 = 

dim I 



(37) 
(38) 

(39) 



7.3 Two-point clustering 

Let us now consider the two-point clustering of the X-ray clusters. 
Our observables will be the angular cross and auto spectra Ct{i, j) 
between clusters binned either in redshift or in photon counts. The 
measured signal is affected by Poissonian shot noise: 



(40) 



where Sij denotes the Kronecker delta. The covariance matrix for 
the angular power spectrum can be written as follows (Hu & Jain 
2004; Cohn 2006; Huterer et al. 2006; Yamamoto et al. 2007): 



S = CoY[Ce{i, j) Cii{m,n)] — Cov^(ij)(mn) 



(41) 



(2£+l)/,ky 



|^C<?(i,m) Ci{j,n) + Ci{i,n) Ci{j,m) 



The covariance matrix is thus block-diagonal, where the number of 
blocks is given by the number of multipoles £ and where the dimen- 
sion of every block is given by the number of distinct pairs which 
can be formed with the adopted number of bins, in redshift or pho- 
ton counts. Assuming a Gaussian likelihood function for the un- 
derlying density field, i.e. for the spherical coefficients aim which 
define the angular power spectrum (aimCier^r) — Ci 5u' 5mm', 
we write the "clustering" Fisher matrix as follows: 



E E 



Cov^ ^ {ij){mn) 



dCi{m, n) 



d9g 



(42) 



8 RESULTS 



8.1 Cluster counts 

We first assume having no information on the cluster redshifts and 
study the constraining power of a number-count experiment. The 
first row in Table 3 summarizes our (optimal) results, obtained with 
about 20 logarithmically-spaced bins in the range 50 < < 10^ 
plus one bin extending to infinity: l-cr errorbars often exceed the 
fiducial values of the parameters themselves and are thus of little 
interest. We only quote errors obtained marginalizing over the en- 
tire parameter set. Note that the conditional errors could be as good 
as A/^L^' = 9 and Aas ~ 10~*. This indicates that covariances 
among parameters are strong (see Fig. 8, upper-left panel). The 
most degenerate pairs are fim-crs, /nl^'-cts, /NL^^'-fim, /3lm- 
(tlm, and olm-cs- It is exactly these covariances that are reduced 
by optimizing the binning scheme in rj, while conditional errors are 
rather insensitive to it. 

Redshift information on the individual clusters must be added 
if we want to constrain the cosmological model and the physics 
of the ICM from a number-count experiment with eROSITA. We 
thus study how the 1-ct errorbars for the various parameters im- 
prove with the accuracy of the redshift estimates. For this calcula- 
tion binning is implemented both in photon counts and in redshift. 
We consider ~ 20 logarithmically-spaced bins in photon counts as 
detailed earlier, while redshift bins have size Az{l + ztin), where 
Az is a parameter we vary and ztin is the median redshift in a bin. 
As soon as some redshift information is available, even as rough as 
with Az ~ 0.1, uncertainties on the parameters shrink by factors 
of ten (e.g. for ag). This applies to the marginal errors, indicat- 
ing that redshift information breaks degeneracies. The conditional 
errors, on the other hand, do not show any dependence on the num- 
ber and size of the redshift bins. In Table 3 we list the uncertain- 
ties that eROSITA should be able to place if a redshift binning 
with Az ~ 0.05 is performed, which approximates what could be 
achieved from photometric estimates ("Counts -I- Photo-z "). 

Spectroscopic-redshift measurements (Az ~ 0.01) would 
result in significantly better constraints (by 30-60 per cent) on 
parameters like n^, h, 51b, Qlm and (tlm with respect to what 
is achievable using photometric redshifts: yet, eROSITA counts 
alone would not improve upon current constraints which have been 
however obtained by combining different probes. Note that the cur- 
rent constraints on the cosmology sector listed in Table 3 refer to 
other probes than galaxy clusters, and the ones on the ICM sec- 
tor are at fixed cosmology. Our results instead have always been 
obtained after marginalizing over both the cosmology and ICM 
sectors simultaneously. Therefore a direct comparison is not com- 
pletely meaningful and the entries "Current Errors" should be sim- 
ply considered as a reference. Note that, if the luminosity-mass 
scaling-relation parameters were exactly known, eROSITA cluster 
counts alone would constrain as and Q^n down to Aag = 0.017 
(2.1 per cent) and Af2m ~ 0.0086 (3.1 per cent; see Table Bl in 
Appendix B). 

Finally, it is worth stressing that, even with some redshift 
information, errors in parameter estimates display strong correla- 
tions. As shown in the upper-right panel of Fig. 8 (for photo-z bins 
and a few tens of bins in photon counts), this is particularly evi- 
dent for the triplet n^, h, and /nl^' within the cosmology sector, 
all the LM parameters among themselves, and for as with all the 
LM parameters (see also Appendix A). 
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Figure 8. Correlation coefficient between pairs of model parameters for tlie diiferent experiments. 
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Figure 9. Choice of i^-nin and fmax in the clustering measurements: l-cr 
errors on the different parameters are shown as a function of £niin (ft fixed 
^max = 96, upper panels) and as a function of ^max (at fixed £niin = 
5, lower panels). The entire eROSITA sample is projected on the same 
sphere and split in 12 bins based on the photon counts rj. 

8.2 Angular clustering 

8.2.1 Ci without redshift information 

Starting from the reference configuration summarized in Table 1, 
we want to optimize the measurement of the cluster 2-point func- 
tion to extract maximal information, for the all-sky survey and with 
multipoles in the range 5 ^ £ ^ 96. We start by assuming no 
knowledge about the redshifts of the individual clusters. The de- 
tected clusters will be binned only in terms of the observed photon 
counts, such that an approximately equal number of clusters are as- 
signed to each bin. Distinguishing in photon counts is fundamental 
to tighten the errorbars both on /nl^' and all the other parameters. 
The gain on /Nf' due to the increment in the number of bins al- 
ready converges after 6 bins: in fact, shot noise is not negligible 
already with a few 10** objects per bin. The best constraints from 
the angular clustering are summarized in Table 3 where we report 
the results obtained using 12 bins ("Angular clustering"): while the 
conditional error on f^l""^ could be as good as d.f^t"^ ~ 24, the 
rather poor constraints on the other parameters increase the corre- 
sponding marginal constraint to ±46 at l-cr. Not surprisingly, the 
non-Gaussianity parameter presents very weak correlations with 
the other ones, once a few bins in photon counts are considered. 
On the other hand, the pairs Us-h and /i-fib are extremely corre- 
lated even with the tightest binning scheme (see Fig. 8, center-left 
panel). Moreover the measurement of the angular power spectrum 
in a single redshift bin cannot break the severe degeneracies which 
plague the ICM parameters, as predictable by simple inspection of 
Fig. 6 (see also Appendix A). 

Note that the minimum multipole which is considered in the 
analysis of the angular power spectrum impacts the estimation of 
the different cosmological parameters (Fig. 9, where the effect of 
^max is also shown). After marginalizing over the whole set of pa- 
rameters (and considering multipoles up to i'max = 96), the uncer- 
tainty on /nl^' improves by a factor of 6 if the minimum multipole 
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is reduced from f ~ 20 to ^ ~ 5. On the other hand, the constraints 
on the other cosmological parameters are insensitive to irain. 

8.2.2 Ci tomography 

How much does the constraining power of Ci measurements im- 
prove if the redshifts of the clusters were available? To answer this 
question we divide the eROSITA sample into redshift bins of size 
Az[l + zun) and, additionally, into 12 bins based on the pho- 
ton counts. In the rest of the paper we will refer to this configura- 
tion as "tomography". By slicing in redshift, the marginal errorbars 
of all parameters drop significantly, even by orders of magnitude. 
Four are the factors which make the tomography more powerful 
than measuring Ci from the full sample: i) the maximum multi- 
pole can be raised to much higher values, reaching £max > 200 
for redshift bins above z ~ 0.8 with fcmax ~ 0.1 h Mpc^^; (ii) 
the Ci signal itself increases for thinner bins due to reduced pro- 
jection effects; (iii) the scale-dependent bias induced by /nl^' is 
more pronounced at high redshift and not averaged out by the pop- 
ulations of lower redshift less-biased objects; (iv) since we always 
consider all cross-spectra among different redshift slices, the num- 
ber of data points scales as n(n + l)/2, n being the number of 
redshift slices. The proper inclusion of all the cross-spectra (wrt to 
the case where only auto-spectra are included) tightens the error 
budget up to 33 per cent (for /nl^') when A2: ^ 0.05 and up to 
40-60 per cent (for f^l'^K n,, h, fib) when Az ~ 0.01. Note that 
even with Az ~ 0.01, there are still more than a thousand objects 
per 2-slice in half of the bins. 

By combining the angular clustering signal from more than 
10 different redshift slices (Az ~ 0.1), it is possible to constrain 
primordial non-Gaussianity of the local type with a marginal er- 
ror < 10 (the corresponding conditional error on the other hand 
is ~ 5 — 6 at 1-a). On the other hand, the gain in the informa- 
tion about fj^l""^ by increasing the number of redshift slices flat- 
tens much faster than almost any other parameter (see also Section 
9.4 and Table 5). In Table 3 we give as a reference the constraints 
obtained with Az ~ 0.05 ("Angular clustering -1- Photo-z "), con- 
sistently with what we adopted for the number-count experiment 
when photometric-redshift estimates are available. For /nl^' we 
find a forecasted marginal error of ~ 10 which is very interesting 
and could potentially rule out entire classes of inflationary mech- 
anisms. With an exposure of 1.6 ks, a tomography experiment is 
more efficient than a "Counts + Photo-z " measurement to con- 
strain the LM parameters, but we have noticed that the situation is 
reversed for higher exposure times (i.e. 3 ks). On the other hand, 
eROSITA tomography alone cannot place constraints on the other 
cosmological and ICM parameters which are by any means com- 
petitive to what is already known nowadays by the combination 
of CMB, BAO and SN data. Note, however, that assuming per- 
fect knowledge of the luminosity-mass relation, the constraints im- 
posed on o"8, and Qm by the eROSITA tomography alone 
would be competitive with the current ones (see Table Bl in Ap- 
pendix B). If viceversa one could assume perfect knowledge of the 
cosmological parameters, tomography alone would largely improve 
our knowledge on the slope and evolution of the LAI relation (see 
in Table 3 the results with Planck priors for an upper limit of the 
uncertainties). 

It is worth mentioning that estimates on /nl'^' based on the Ce 
tomography are nearly perfectly uncorrelated from any other model 
parameter (Fig. 8, center-right panel) although strong correlations 
are still present in the complementary parameter space. 

Fig. 10 allows us to stress the importance of high-redshift ob- 
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Figure 10. Marginal error on /J^'^"' and total Figure of Merit from different 
individual redshift bins. The bin width is A2(l + z) with Az = 0.05 for 
a total of 16 bins up to z ~ 1.5. Magenta circles indicate the constraints 
on /^L^' obtained from every individual redshift slice, the horizontal ma- 
genta line is the global constraint on obtained considering all the 
bins together, including cross-correlations among redshift bins. The Figure 
Of Merit of Eq. (34) is color coded in shades of gray for the individual z- 
bins, while the total FoM of the tomographic measurement is indicated in 
the upper-right box. 

jects: there we show how the total figure of merit (as defined in 
Eq. (34) and which takes into account the errorbars and corre- 
lations among the whole set of 10 parameters) and the marginal 
constraints on /Nf' depend on the individual 2— bin, within our 
binning scheme; the magenta horizontal line refers to the marginal 
constraint on /nl^' given by the combination of all redshift slices, 
including cross-correlation among bins; also the number of objects 
in each redshift bin is indicated. The final constraints are the re- 
sults of at least three factors: the thickness and the redshift of the 
z-bin, and the number of objects in each bin. Considering only the 
slice with 0.76 < z < 0.85 (which contains ~ 2800 clusters) 
gives already A/^l^' = 51. In general, the sample of clusters with 
z > 2med better constrains the cosmological parameters (by a fac- 
tor of a few) wrt the identically sized sample with z < Zmcd ■ 

8.3 Combining counts and clustering 

We want now to combine our results on the cluster counts and on 
the angular power spectrum to get tighter constraints on the model 
parameters. This requires accounting for the cross covariance be- 
tween the two probes, which is proportional to the angular bispec- 
trum of the clusters (see e.g. Takada & Bridle 2007). For Gaussian 
initial conditions the bispectrum vanishes on large scales. Cross 
covariances remain small also for the weakly non-Gaussian mod- 
els we consider, in spite of the large cluster bias factors. This is 
because we only consider very large scales and rather heavily pro- 
jected data. For this reason we ignore cross covariances between the 
cluster number counts and the measurement of the angular power 
spectrum. As a result, the total Fisher matrix is obtained by sum- 
ming up the single Fisher matrices of the individual experiments. 
The results obtained considering the joint datasets are summarized 
in Table 3, in Fig. 11, where error contours for a selection of pa- 



rameter pairs are shown, and in Fig. 8 where we quantify the cor- 
relations among parameter pairs for all the different probes. Once 
again, we distinguish the cases in which redshift estimates for the 
individual clusters are available or not. For all parameters, measur- 
ing photometric redshifts (with Az ~ 0.05) shrinks the marginal 
errorbars by a factor of 4 or 5. Adding information on the num- 
ber counts to tomographic Ce measurements has little impact on 
/nl''' but significantly shrinks the confidence interval of the other 
cosmological and ICM parameters. On the other hand, "Counts -l- 
Photo-z " alone gives better results than "Counts -I- Angular clus- 
tering" without photo-z for all parameters but /nl^' arid ris- Prop- 
erly taking into account the uncertainties on the LM parameters 
strongly degrades the cosmological constraints which could be ide- 
ally achieved if the cluster observable-mass mapping was perfectly 
known (e.g. "Counts -I- Angular clustering -I- Photo-z" vs "Counts 
-I- Angular clustering -I- Photo-z -I- LM fixed", and Table B 1). 

8.4 Adding external priors 

X-ray cluster counts and power spectra (with and without photo- 
metric redshifts) cannot simultaneously determine all the model 
parameters (see Table 3 and Fig. 11): in particular the triplet Us, 
h, and f2b are left unconstrained by eROSITA data. We there- 
fore complement eROSITA with strong standard priors on the 
Hubble parameter (from the Hubble Key Project Ah — 0.08, 
Freedman et al. 2001) and on the mean baryon density of the 
universe (from primordial nucleosynthesis, A(Slb h^) = 0.002, 
Kirkman et al. 2003), as commonly done in X-ray studies (e.g. 
Mantz et al. 2010b). The sensitivity of our results to these priors 
is shown in Table 3 ("-I- Priors "). The parameter that benefits the 
most is the spectral index, while all the other cosmological and 
ICM parameters are rather insensitive to them (at least for the joint 
constraints with photo-z). While constraints on h and fib are domi- 
nated by the imposed priors for the abundance experiments and the 
clustering without redshifts, eROSITA's tomography (and thus the 
combinations of the experiments) significantly contributes to fur- 
ther shrink those errors. 

We focus now on the CMB analysis performed by the Planck 
satellite whose results will be available by the time the eROSITA 
all-sky survey will be completed. We consider the Planck Fisher 
matrix for a measurement of the power spectrum of temper- 
ature anisotropies calculated by the Dark Energy Task Force 
(Albrecht et al. 2009), kindly made available to us by Dragan 
Huterer and Wayne Hu. This constrains all cosmological param- 
eters but /nl^' (we do not consider the CMB bispectrum in this 
work). The results obtained using the Planck forecast as a prior for 
our eROSITA analysis are shown in Table 3: the constraints on 
all parameters are strongly improved, in particular for Ua, h, flh 
(where the improvement reaches a factor of 30!). Yet, eROSITA 
significantly contributes to the determination of all parameters. 
Cosmological constraints on ag (which is not a "natural" choice 
to fit CMB data) and r^m will be known at ~ 1.5 — 2 per cent 
accuracy, while constraints such as Af^£^^ ~ 7 will be compara- 
ble to the results of the CMB three-point statistics (but on different 
spatial scales). Similarly, all the ICM parameters will be more ac- 
curately determined than they are nowadays (see error ellipses in 
Figs. Al, A2, and A3). Note that in terms of LM parameters only, 
adding Planck priors to the cosmological sector ("Counts -I- Angu- 
lar clustering -I- Planck") is almost equivalent to measuring pho- 
tometric redshifts for the whole cluster sample ("Counts -I- Angu- 
lar clustering -I- photo-z"). A comparison between the constraining 
capabilities of Planck and eROSITA is given in Table B2 in Ap- 
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Table 3. Marginal 1 — cr errors for the cosmological and the LM parameters obtained using the Fisher-matiix formalism. The label "+Priors" indicates results 
obtained adopting external priors on h and as detailed in Section 8.4. The shaded row marks the best constraints obtained using eROSITA data only. We refer 
to the following survey: Toxp = 1.6ks, rj^in = 50, M^in = 5 x IQ^^h'^MQ, /^ky = 0.658, Az = 0.05. 



eROSITA data 


EoM 




Acts 




Ans 


Ah 


Afib 






A/9lm 


A(7LM 


Counts 


1.0 


~ 9 X 10^ 


~ 1.6 


~ .5 


^ 4 


^ 4 


~ .3 


~ 1.8 


~ 7 


~ 9 


~ 3 


Coiints + Priors 


4.6 


~ 8 X 10^ 


~ 1.5 


~ .4 


~ 2 


.080 


.0113 


~ 1.7 


~ 7 


~ 7 


~ 2 


Counts + Photo- 2 


10.7 


423 


.113 


.0191 


.559 


.558 


.0649 


.20 


.75 


~ 1 


.277 


Counts + Photo- 2 + Priors 


12.8 


360 


.100 


.0188 


.205 


.078 


.0110 


.16 


.73 


~ 1 


.202 


Angular clustering 


7.1 


46 


.257 


.0817 


.845 


~ 1 


.0974 


.47 


~ 1 


~ 3 


~ 1 


Angular clustering + Priors 


9.1 


42 


.226 


.0693 


.256 


.068 


.0100 


.36 


~ 1 


~ 3 


~ 1 


Angular clustering -f Photo-2 
























Angular clustering + Photo-2 + Priors 


13.2 


9.8 


.095 


.0207 


.076 


.028 


.0033 


.08 


.43 


~ 1 


.242 


Counts + Angular clustering 


10.6 


42 


.180 


.0582 


.530 


.967 


.0736 


.24 


~ 1 


~ 2 


.621 


Counts + Angular clustering + Priors 


12.5 


37 


.169 


.0531 


.154 


.064 


.0089 


.22 


~ 1 


~ 2 


.557 


Counts + Angular clustering + Photo-2 


16.3 


8.8 


.036 


.0118 


.088 


.153 


.0114 


.05 


.20 


.397 


.117 


Counts + Angular clustering -t- Photo-2 + Priors 


17.2 


8.2 


.036 


.0111 


.033 


.027 


.0030 


.05 


.20 


.384 


.114 


Counts + Angular clustering -t- LM fixed 




36 


.016 


.0099 


.172 


.461 


.0464 










Counts + Angular clustering -t- Photo-2 + LM fixed 




8.4 


.003 


.0029 


.055 


.110 


.0092 










Counts + Angular clustering -f Planck 


19.7 


26 


.022 


.0065 


.004 


.005 


.0005 


.04 


.16 


.348 


.137 


Counts + Angular clustering -f Photo-2 + Planck 


22.1 


6.9 


.014 


.0039 


.003 


.003 


.0003 


.02 


.07 


.173 


.045 


Current Errors'^ 




[-10,+74] 


.024 


.0061 


.012 


.014 


.0016 


.14 


.42 


.085 


.039 


Planck Errors 






.024 


.0071 


.004 


.006 


.0006 











^ WMAP7+BAO+SN for the Cosmology sector independent of the ICM sector, Komatsu et al. (201 1 ); results by Vikhlinin et al. (2009a) for the LM parameters, 
at fixed Cosmology. 95.4 per cent CI for /^l^' 

^ Error estimates for future power- spectrum measurements with Planck: marginalization solely over the cosmology sector, excluding /JJl^^ 
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Figure 11. Joint I-ct error ellipses for a selection of parameter pairs obtained by marginalizing over all the other model parameters (no priors assumed). 
Individual experiments and their combinations for the all-sky survey (described in Table 1) are indicated with different colors. Contours for the number-count 
experiment with no redshift information vastly exceed the area in the plot. The complete set of eiTor ellipses is given in Appendix A, with priors oh h and H^,. 
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pendix B, where the analysis is repeated exclusively for the cosmol- 
ogy sector and assuming Gaussian initial conditions (5-parameter 
fit). Note that, if we could perfectly characterize the X-ray cluster 
scaling-relations, eROSITA would outperform Planck in the deter- 
mination of (78 while, even in this idealistic scenario, Planck would 
still do significantly better for ris, h, and fib- 



9 DISCUSSION 

9.1 Wide or deep surveys? 

Is there any advantage in undertaking a deeper survey than the all- 
sky one analyzed so far? We consider here an hypothetical survey 
with average exposure time of 7.5 ks. Accordingly, we assume a 
sky coverage of 6,000 deg^ (with imin ~ 9), which gives approx- 
imately the same amount of observing time used for the all-sky 
survey. Although the higher exposure time enables the detection of 
fainter clusters at high redshift, the limited sky coverage signifi- 
cantly reduces the total number of sources. With a detection limit 
of 50 photons and a minimum mass of 5 x 10^^ K'^Mq, the deeper 
and less extended survey should identify 6.85 x 10* clusters with 
a median redshift of 0.56. 

These changes impact the cosmological constraints in differ- 
ent ways. For a number-count experiment, the deeper survey better 
constrains some parameters such as /nl*'> «s, and 7lm (15-30 per 
cent improvements) but is not optimal overall. On the other hand, 
clustering studies clearly benefit from the all-sky coverage (param- 
eter constraints improve by many tens per cent), especially if one 
wants to measure /nl*' through the scale-dependent bias which is 
only evident on large projected angular scales. In Table 4, we re- 
port the total number of detected clusters and some constraints on 
the model parameters that can be obtained from the joint analysis 
of cluster counts and Ci tomography (both with redshifts of photo- 
metric quality) for a series of survey strategies, where the covered 
sky fraction and the average exposure time are changed with re- 
spect to our reference choice. For example, the 3 ks survey over 
27 '000 deg^ could be realized if, instead of performing three years 
of pointed observations after the first four years of the eROSITA 
mission, the all-sky survey would be continued for the remaining 
lifetime of the satellite. This, of course, would have the drawback 
that the planned pointed cluster follow-up observations could not 
be performed. The mild improvement in the resulting constraints 
might not justify this possible extension of the all-sky survey. 

9.2 Does 1/2+1/2 equal 1? 

It has been agreed that eROSITA data will be equally split be- 
tween the Russian and the German consortia: each of them will 
own an equal fraction of the sky. We want to investigate what would 
be the impact of performing separate analyses for the two half-sky 
surveys. Two effects have to be considered: halving the number of 
objects and the loss of the largest angular spatial separations where 
the effect of primordial non-Gaussianity is larger. As an example, 
we report how passing from /gky = 0.658 to /sky = 0.33 (with 
tmin ~ 7) impacts the constraints on /nl*' and as. We find that 
the joint marginal errors for the all-sky survey A/^l''' ~ 8.8 and 
Acts = 0.036 degrade to A/i°£='' = 18.9 and Aag = 0.062. 
Combining the results from the two halves of the sky a posteriori 
would then give: A/^l"' = 13.4 and Acts = 0.044 (where we 
have optimistically assimied that the two halves of the sky are sta- 
tistically independent). 



eROSITA Data 


FoM 


A/iff"' 


Acts 


Counts 


1.1 


~ 9 X 


~ 1.6 


Counts + Photo- z 


10.7 


423 


.113 


Counts + Spectio-2 


11.2 
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.095 


Angular clustering 


7.1 


46 


.257 


Angular clustering + Photo- z 


12.0 


10.1 


.097 


Angular clustering + Spectro-z 


14.5 


7.8 


.059 


Counts + Angular clustering 


10.6 


42 


.180 


Counts + Angular clustering + Photo-z 


16.3 


8.8 


.036 


Counts + Angular clustering + Spectro-z 


17.7 


7.0 


.024 



Table 5. Forecasted total figure of merit, and 1-<t marginal errors on /JJ'l'*' 
and ITS obtained with the eROSITA cluster sample. First, no redshift in- 
formation is assumed. Then the data are sUced in redshift bins of width 
lS.z{l + z), with Az = 0.05 and 0.01 (mimicking photometric vs. spectro- 
scopic redshift estimates). No priors assumed. 



9.3 The case with no clusters beyond z ^1 and 
"the magnificent 1000" 

Many parameters of interest exhibit stronger effects on the clus- 
ter number coimts and Cg. at higher redshifts than at lower red- 
shifts: this is particularly pronounced for /nl°^'. but also evident 
for some of the ICM parameters. Measurements at high redshifts 
may be difficult or even out of reach for at least two reasons: i) the 
higher the redshift, the more problematic the AGN contamination 
might be; ii) the higher the redshift, the lower the chances of suc- 
cessfully measuring a (photometric or spectroscopic) redshift are. 
We briefly compare here the results discussed so far with the pes- 
simistic scenario where no clusters beyond redshift z ^ 1 can be 
included in the analysis. The total number of objects which would 
be lost beyond redshift 1 is about 1200 (for Texp = 1.6 ks; three 
times higher with a double exposure time). With redshift bins of 
width 0.05(1 -I- z), the constraints on /nl''' would deteriorate to 
A/nl''' = 517 and Af^t^ = 11.1 for a number-count and a to- 
mography experiment, respectively (see also Table 4 for the joint 
constraint). Similarly, the loss of accuracy in measuring all the 
other parameters would be more pronounced for number counts, 
with relative deteriorations of about 10-20 per cent. 

Let us now imagine to only use the most massive 1000 clus- 
ters at 2: ^ 1 for our analysis ("the magnificent 1000"). This sample 
includes objects with masses above ~ 2.2 x W^^h~^ Mq, possi- 
bly similar to what Sunyaev-Zel'dovich surveys might yield in the 
future. The joint analysis of number counts and Ce tomography 
with this subsample would give a marginal error of A/^l^' ~ 41 
(with photometric-redshift information). The attractive conditional 
error of Af^^'^^ ~ 12 is in practice uninteresting because plagued 
by huge systematics and parameter degeneracies: with the "mag- 
nificent 1000" only the exponential cut off of the mass function 
is probed, where the models are more uncertain and where ev- 
ery parameter is in practice degenerate with another Note that 
for the analysis with the whole eROSITA sample, the compari- 
son between conditional and marginal error is much less dramatic, 
reading Af^l""^ ~ 4.4 vs A/Jj^"' - 8.8 (see Table 3). In any 
case, the marginal constraints on all the other parameters would 
be much less appealing (Table 4). Combining the analysis of the 
magnificent 1000 with Planck priors would still give A/^l''' ~ 26 
and Aos = 0.024 (corresponding to a relative error of 2.9 per 
cent), a figure which is destined to further improve with multi- 
wavelength follow-ups of the cluster sample and consequently bet- 
ter mass proxy measurements. 
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8.6 
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1.6 


4.66 X 10'' 


18.9 


6.5 


4.6 


(111 + 1/2) eROSITA 


27'000 


1.6 


9.32 X lO-* 


13.4 


7.3 


5.1 


ziil eROSITA 
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1.6 


9.19 X 10* 


9.5 


7.8 


5.2 


"The magnificent 1000" 


27'000 


1.6 


~1000 


41 
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Table 4. Forecasted performance of different survey strategies (no priors assumed). The reference case (labelled eROSITA) refers to the all-sky survey 
extensively discussed in Table 3. In all cases we consider a detection limit of 50 photons and a minimum cluster mass of 5 X 10^^ h~^MQ. Results are given 
for the joint analysis of number counts and angular clustering with photometrically determined cluster redshifts, i.e. the data are split in redshift bins of width 
0.05(1 + z). We use here separate figures of merit to distinguish the cosmological parameters from those of the LM relation. 



9.4 Spectroscopic redshifts 

Our reference forecasts assume that photometric redshifts of the in- 
dividual clusters will become available. However, given that spec- 
troscopic follow-ups of the eROSITA survey are being proposed, 
in Table 5 we present new results that refer to the optimistic sce- 
nario where spectroscopic redshifts for all the detected clusters 
will be available. In this case we use redshift bins of thickness 
0.01(1 + z). For all model parameters improvements are much 
more pronounced for the Ce tomography rather than in a count 
experiment. Marginal constraints on all the cosmological param- 
eters but /nl'^' ™d exhibit an improvement of more than 50 
per cent when passing from photometric- to spectroscopic-redshift 
quality with a tomography measurement alone; all parameters but 
fNL^^ get constrained 30-40 per cent tighter when joint constraints 
with spectroscopic redshifts are considered wrt to the photo-z case. 
Overall, the total Figure of Merit would improve substantially if 
spectroscopic redshifts were available. 

9.5 The mass cut 

So far we have always applied a cut in the -q-z plane which cor- 
responds to a cluster mass of A/ = 5 X 10^^ h-^AlQ in the 
fiducial model. Moving this threshold significantly changes the se- 
lected cluster population both in terms of their sheer numbers and 
redshift distribution (see Table 2 for details). Of course there are 
more objects when the mass threshold is lowered and this for- 
mally produces better parameter constraints from counting exper- 
iments. On the other hand, a lower threshold reduces the effective 
bias of the cluster population (and thus the non-Gaussian correc- 
tion) and moves also down the median redshift (and thus fmax 
for the clustering measurement): for the clustering experiments, 
trends are different for different parameters. The joint constraints 
for /nl^' are rather insensitive to the mass threshold, while the pa- 
rameters which are most affected by the exact value of the mass 
cut are erg and 7lm (their constraints improve by 30 per cent when 
Mcut = 10" H-^Mq). 

A note of a caution is in order here. In a meaningful Fisher- 
matrix study the sample definition (and thus the mass cut) must not 
change with the model parameters. We identify the locus in the r)-z 
which corresponds to a cluster with the threshold mass in the fidu- 
cial model and never modify it when we compute the derivatives 
with respect to the parameters. This mimics what would be done in 
practice to select high-mass objects out of an observational sample. 
The sample selection is done before computing the statistical quan- 
tities. For this reason it would be incorrect to impose a cut based 
on the actual mass or X-ray luminosity of the clusters as sometimes 
it has been done in the literature (unless mass or luminosity are a 



direct observable, which is not the case for eROSITA). We have 
checked that adopting this procedure leads to unrealistically opti- 
mistic forecasts. 



10 OTHER NON-GAUSSIAN MODELS 

Primordial non-Gaussianity of the local type is only one (the 
simplest) of the countless ways to perturb the primordial grav- 
itational potential around the Gaussian assumption (e.g Chen 
2010). Here we present the error budget that eROSITA should 
achieve shall the data be fitted with two other models/templates 
for early non-Gaussianity, the so-called orthogonal and equi- 
lateral types. These models differ from the local one in the 
sense that the bispectrum of the Bardeen's potential peaks at 
different triangle configurations and are characterized by a sin- 
gle parameter, /nl'^° or /^l" ' which quantifies the ampli- 
tude of such bispectrum (Senatore, Smith & Zaldarriaga 2010; 
Creminelli et al. 2006). While the local shape is expected in 
all models where non-linearities develop outside the horizon 
(e.g. in multi-field inflation), equilateral non-Gaussianity is pro- 
duced in single-field inflationary models with a non-minimal La- 
grangian where the mode coupling is created by non-canonical 
kinetic terms or higher-derivative operators (e.g. in ghost infla- 
tion or DBI^* inflation). Finally, the orthogonal template repro- 
duces (among the others) non-Gaussian features arising from as- 
suming non-standard vacuum choices. These three shapes are "or- 
thogonal" to each other (Babich, Creminelli & Zaldarriaga 2004; 
Senatore, Smith & Zaldarriaga 2010) and linear combinations of 
them can arise in physically-motivated models of the early uni- 
verse. 

In terms of late-time observables, the striking difference be- 
tween these models lies in the scale-dependence of the bias for 
DM haloes: at large separations, the bias scales as ' 
and f^i}^°k^^ for the local and the orthogonal types respec- 
tively, while it is asymptotically constant (but, for fixed halo 
mass, numerically different than in the Gaussian case) for the 
equilateral model. The different shapes of the primordial bispec- 
tra imply differences in the skewness of the linear density field 
and, consequently, different modifications to the Gaussian halo 
mass function too, via Eqs. (8) and (9). We adopt the formu- 
lae for the skewness from Taruya, Koyama & Matsubara (2008) 
and Schmidt & Kamionkowski (2010), and the improved formulae 
for the non-Gaussian halo bias by Desjacques, Jeong & Schmidt 
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Table 6. Marginal l-cr constraints enabled by eROSITA for different models of primordial non-Gaussianity. Assumptions and the survey strategy are as in 
Table 3 with no external priors unless explicitly stated. 



(2011): these have shown better agreement with N-body simula- 
tions than previous analytical derivations. However, these expres- 
sions are only valid asymptotically as fc — )■ and cannot be used at 
much larger k where the halo bias is basically constant and the pro- 
vided formulae do not reproduce such behavior. In order to solve 
this problem, we have manually imposed this asymptotic behav- 
ior for large k by removing the scale dependencies due to the lin- 
ear transfer function T{k) at the wavenumbers where the factor 
1/T(k) starts to dominate. 

As detailed in Table 6, given the less prominent scale- 
dependence of the halo bias, the constraints on ./nl*^" and /^l"'' 
are weaker than for the local case; the constraints on the other cos- 
mological parameters and those for the LM relation also exhibit 
a non-negligible modification. Counts and tomography joint con- 
straints read 0.053 and 0.049 for erg in the orthogonal and equi- 
lateral models, respectively, instead of 0.036 with the local tem- 
plate. LM parameters worsen by 10-20 per cent and 20-40 per cent 
with the orthogonal and equilateral models, respectively, while the 
triplet Us, h, 0,^ improve in both cases by a few per cents. 



11 ADDITIONAL REIMARKS 

In this Section we critically discuss some of the assumptions and 
the methods we have been using to derive our main results. 

11.1 Theoretical models for the halo mass function 

Since the parameters of the mass function are determined from N- 
body simulations, they are known with some intrinsic uncertainty 
that we have ignored in our main study. One possible way of gaug- 
ing this uncertainty is to introduce in the Fisher-matrix analysis a 
set of nuisance parameters (as many as the mass-function parame- 
ters) with their own prior distribution (which embodies the covari- 
ance matrix from the fits to the simulation) and marginalize over 
them (see Wu, Zentner & Wechsler 2010; Cunha & Evrard 2010, 
for dark-energy surveys). We have implemented this scheme and 
repeated our calculations using the covariance matrices for the pa- 
rameters Aj^^, o™^, 6™'' and c™^ kindly made available to 
us by Jeremy Tinker. Unfortunately, owing to the procedure used 
to fit the simulations, a covariance matrix including also the 3 evo- 
lutionary parameters A^^^jO^^^ and a is not available. Any- 
way our analysis should at least give an order-of-magnitude esti- 
mate of the effect. The errors provided by Jeremy Tinker on his 
z=0 mass function parameters at fixed cosmology translate in an 
overall uncertainty of 1 per cent in the shape of the mass func- 
tion: it is the variation in Cosmology which is responsible for the 



5 per cent accuracy level quoted in Tinker et al. (2008). In order to 
approximately take into account possible larger uncertainties, we 
have rescaled the covariance matrix so that to obtain a marginal 
uncertainty of about 10 per cent in the abundance of objects with 
M ~ 10^^ h~^MQ at 2: = before using it as a prior on the mass- 
function parameters. We find that, within this setup, the constraints 
on the parameters both of the cosmology and ICM sector are de- 
graded only by up to 10 per cent when 4 additional mass func- 
tion parameters are included in the fitting analysis and marginal- 
ized over. These constraints are destined to further degrade when 
adding the 3 evolutionary parameters. IVIoreover, this picture does 
not take into account possible systematic errors (up to 20 per cent) 
that baryonic physics seems to induce in the halo mass function 
(Stanek, Rudd & Evrard 2009; Cui et al. 201 1). 

Imperfect knowledge on the large-scale bias of the DM haloes 
should have a smaller effect and will not be discussed here (see 
Wu, Zentner & Wechsler 2010). 

11.2 IVIass-luminosity relation 

Our forecasts rely on a series of assumptions for converting DIVI 
halo masses into X-ray photon counts. 

1) We have assumed a power-law LM relation. While there is, 
as of now, no very compelling evidence that this assumption breaks 
down (see e.g. Sun et al. 2009; Eckmiller, Hudson & Reiprich 
2011), the scatter has been shown to increase for low-mass clus- 
ters (Eckmiller, Hudson & Reiprich 201 1). For this reason we have 
introduced a minimum mass threshold and varied it within a rea- 
sonable range to test the stability of our results. 

2) We have assumed that the scatter around the power-law 
relation is lognormal, but no consensus has been reached regard- 
ing the amplitude and nature of this noise component (see e.g. 
Vikhlinin et al. 2009a; Pratt etal. 2009; Stanek etal. 2010). As 
shown by Shaw, Holder & Dudley (2010), an unknown amount of 
asymmetry in the scatter of the observable-mass relation can affect 
cosmological constraints in the same manner as additional uncer- 
tainties in the fiducial value of the scatter ctlm itself. 

3) We have used the current limited knowledge on the LM 
relation to fix the parameters that define our fiducial model. For 
instance, we used a scatter in the LM relation of 40 per cent 
(Vikhlinin et al. 2009a). However, more recent observational scal- 
ing relations for core-excised clusters show a reduced dispersion, 
although spanning somehow smaller redshift ranges (IVIantz et al. 
2010a). We have repeated our analysis adopting the evolving 
scaling-relations measured by IVIantz et al. (2010a), converted in 
the (0.5-2.0) keV band: we find that the number of clusters detected 
by eROSITA increases by more than 30 per cent, with a total pop- 
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ulation of 1.24 x 10^ objects (more massive than 5 x 10^^ H'^Mq) 
and a slightly lower median redshift, ^median ~ 0.31. Conse- 
quently, the joint constraints on all the cosmological parameters 
are improved by 10-50 per cent, with the exception of /nl''' ^rid 
fini which show a smaller improvement. On the other hand, the 
marginal constraints on the LM parameters themselves (but not on 
the LM scatter) are tighter by a factor of a few. The choice of the 
fiducial prescription for the LM relation changes more the predic- 
tions for the cluster counts than for the clustering measurements, 
and /nl^' is the parameter which exhibits the smallest dependence 
on it. In any case, also in the light of our (possibly conservative) 
predictions, one of the major results of this paper is that eROSITA 
and Planck results combined together will be able to constrain the 
slope and the evolution factor of the LM relation at per cent level 
(see Table 3). 

4) We assumed that the slope and the scatter of the LM re- 
lation do not depend on redshift. This is basically Occam's razor 
applied to current data but might not be adequate for eROSITA 
accuracy. 

5) We always assumed no prior knowledge for the measure- 
ments of the ICM parameters which, in our study, are fully de- 
termined by eROSITA data. Adding priors in fact would require 
assuming a covariance matrix for the parameters. Results adopting 
perfect prior knowledge of such parameters are given as a reference 
in Table Bl. 



11.3 Temperature-mass relation 

Fig. 1 and the discussion of Figs. 5 and 6 not only show that the 
photon counts depend very weakly on the X-ray cluster temperature 
but also that the exact TM relation for galaxy clusters has a little ef- 
fect on our observables. We have checked anyway what happens if 
4 additional parameters regulating this scaling relation are included 
in the analysis of Table 3 (10-1-4 parameter fit). We found that this 
leads to a deterioration of the constraints in the cosmology sector 
(by an amount ranging from a few per cent to a factor of a few for 
/nl''') as the new parameters are poorly constrained by the data. It 
is worthwhile mentioning that the deterioration on the cosmologi- 
cal constraints (6 parameters) due to 4 additional temperature-mass 
scaling-relation parameters is up to three times less severe than the 
effect of 4 additional parameters regulating the LM relation and its 
redshift evolution (see Appendix B and Table B3). 



11.4 Limber approximation 

The Limber approximation is often used to forecast the constrain- 
ing power of the angular power spectrum for tracers of the cosmic 
large-scale structure. It is well known, however, that this approxi- 
mation, for Gaussian initial conditions, becomes less and less accu- 
rate when projections are taken within thinner and thinner redshift 
bins (Simon 2007; Loverde & Afshordi 2008). We find that this ef- 
fect becomes even more important in the presence of primordial 
non-Gaussianity of the local type or whenever a scale-dependent 
DM halo bias changes the shape of the angular power spectrum on 
large scales. With /nl''' ~ 100, the Limber approximation devi- 
ates from the exact calculation by up to 20 per cent for ^ ~ 5. 
Remarkably, the sign of the discrepancy reverts depending on the 
actual value of /nl^' (see also Giannantonio et al. 201 1). Therefore 
the inaccuracy of the Limber approximation is particularly damag- 
ing to compute the Fisher matrix for a clustering experiment where 
derivatives of Ci with respect to the individual parameters must 



be taken. We have noticed that using the Limber approximation all 
over the considered multipole range to predict the performance of a 
eROSITA Ct study can underestimate the marginal error on /Jjl^' 
by a factor of 4. 

11.5 Tlie mildly non-linear regime 

As discussed in Section 8.2, the current accuracy to which we 
model non-linearities both in the matter-density field and in the 
cluster bias prevents us to trust predictions including information 
from scales with k > fcmax — 0.1 h Mpc~^. We have formally 
tested that if one could (with more robust models for the mildly 
non-linear regime) increase this limit to fcmax ~ 0.3 h Mpc~^, 
the joint eROSITA constraints (with photometric redshifts) might 
tighten by 20-50 per cent for all parameters but /nl^'- 

11.6 Other issues 

In our analysis, the cosmology and ICM sectors are not fully inde- 
pendent. In fact, Eq. (18) introduces a correlation between the LM- 
evolution parameter and the cosmology sector through the term 
7LM \nE{z) . On the other hand, in Section 1 1 .2, we have presented 
forecasts based on the alternative LM relation by Mantz et al. 
(2010a), where the evolution scales as 7lm logio(l + z) and does 
not couple with the cosmological parameters. Also Jib should cor- 
relate with the LM relation since the cosmic baryon fraction relates 
to the total gas mass of a cluster and thus to its observed X-ray 
brightness. We postpone to future work the choice of a proper pa- 
rameterization for this effect. 

Throughout our study, we have assumed a detection limit for 
the galaxy clusters of 50 photon counts, assuming the total flux 
gets detected. We defer to future work the exploration of different 
selection criteria, e.g. based on the signal-to-noise ratio. Moreover, 
we will address the impact of changing the photon-count limit and 
compare our findings with the output of the eROSITA sub-sample 
with strong observed fluxes (say above 1000 photon counts) and 
thus better mass proxies than Lx ■ 

Finally, we have assumed a constant exposure time Texp 
and a uniform hydrogen column density Nh, instead of tak- 
ing into account the possible variations of these quantities across 
the survey area. A more detailed treatment of the distribution of 
Texp and Nh may introduce corrections to our forecasts: for ex- 
ample, within Sunyaev-Zel'dovich cluster count measurements, 
Khedekar, Majumdar & Das (2010) have shown that the combina- 
tion of a wide with a deep survey, keeping fixed the total exposure 
time, gives better constraints than a longer-exposure wide survey 
only. Further study is required to assess whether this is the case 
also for X-ray cluster count and clustering experiments, and when 
fluctuations of Texp and Nh along different lines of sight are not 
averaged over large patches of the sky. 



12 CONCLUSIONS AND SUMMARY 

We have forecasted the accuracy with which the X-ray telescope 
eROSITA will constrain cosmological parameters based on the 
number counts and spatial distribution of galaxy clusters. Our ref- 
erence case considers 10 model parameters 6 of which charac- 
terize the cosmological model (vanilla, flat, ACDM model) while 
the remaining 4 describe the luminosity-mass relation for galaxy 
clusters. Special attention is dedicated to the primordial non- 
Gaussianity parameter /nl^'. 
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Galaxy clusters are sorted in terms of the photon counts that 
will be detected by the eROSITA telescope. We convert the masses 
of DM haloes (for which we can predict abundance and clustering 
properties as a function of the cosmological parameters) into this 
observable by taking into account observationally-motivated scal- 
ing relations, the properties of the X-ray detector, and the integra- 
tion time of the observations. We find that, in an all-sky survey with 
a typical exposure time of 1.6 ks, eROSITA will observe 9.3 x 10^ 
galaxy clusters (more massive than 5 x 10^^ Mq) above the de- 
tection limit of 50 photons in the (0.5-2.0) keV band. Their redshift 
distribution will be broad, extending up to z ~ 1.5 with a median 
of 0.35. 

Our forecast is based on measuring the abundance of galaxy 
clusters as a function of the X-ray photon counts and the corre- 
sponding angular power spectra, via a Fisher matrix approach. We 
combine the two experiments without assuming any prior knowl- 
edge of the X-ray scaling-relation parameters (the so-called self- 
calibration technique). We distinguish two cases based on whether 
redshift information on the individual clusters is available or not. 
Finally, we integrate eROSITA data with priors on the Hubble 
parameter and on the mean baryon fraction, and with future con- 
straints from the Planck satellite; we analyze different survey strate- 
gies, study the impact of model uncertainties, and explore how con- 
straints degrade when extra nuisance parameters are considered. 
Our main findings are reported in Tables 3, 4, 5, and 6. They can be 
summarized as follows. 

• Despite the unprecedented size of the eROSITA all-sky sam- 
ple, without any knowledge of the redshifts of the individual clus- 
ters, it is not possible to simultaneously improve currents con- 
straints on cosmology (with the exception of /nl^') and on the LM 
scaling relation. Redshift information is vital to break the strong de- 
generacies among the parameter estimates. The availability of pho- 
tometric redshifts, with an accuracy better than 0.05(1 + 2), already 
shrinks the constraints on all parameters by a factor of many times 
(see Table 3). 

• Binning the data in X-ray photon counts (77) already removes 
some degeneracies but, in general, binning in redshift is more effi- 
cient in tightening the errorbars on the model parameters, although 
much more expensive in terms of observing time. Anyway, binning 
in T] should be always adopted to achieve the optimal constraints 
independently of the accuracy to which the cluster redshifts are 
known. 

• A tomographic study of the angular power spectrum with at 
least 10 redshift slices (i.e. with < 0.1(1 + 2), where the data 
are binned also in photon counts) gives exceptional l-o marginal 
errorbars on the non-linearity parameter of A/^l^' < 10. This 
also provides tighter constraints than any number-count experiment 
in (photometric-) redshift bins for all cosmological parameters with 
the exception of fim. The combination of 1- and 2-point statistics in 
tomographic slices based on photometric redshifts is optimal, with 
A/i!j°£='' ~ 9, Acrg ~ 0.036 (4.4 per cent) and AS^Jm ==; 0.012 (4.2 
per cent) (see Table 3). 

• Also for the LM scaling relation, Ci tomography gives tighter 
constraints than cluster number counts, although this result is re- 
versed with a deeper exposure (Tcxp = 3 ks). The joint analysis of 
the two eROSITA experiments improves the current errorbars on 
the slope of the LM relation by a factor of 3 (see Table 3). 

• Measuring spectroscopic redshifts for the individual clusters 
would further tighten the marginal error on each parameter by an 
additional ~ 30 — 40 per cent wrt to the photo-2 case. This ap- 



plies to all parameters but /nl^', which exhibits somewhat smaller 
improvements (see Table 5). 

• Yet, even the optimal combination of eROSITA measure- 
ments with redshift information only loosely constrains parameters 
such as Us, h, and Sib- Standard priors on h and Sib are necessary 
to reduce uncertainties for this subset of the parameter space (see 
Section 8.4). 

• If we could perfectly characterize the LM scaling relation 
with prior data, eROSITA (with photometric redshifts) would con- 
strain the vanilla ACDM in line with the future CMB mission 
Planck, with unrivaled constraints on the amplitude of the linear 
DM power spectrum erg (to 0.4 per cent accuracy, see Tables Bl 
and B2). In fact, while accounting for primordial non-Gaussianity 
does not sensibly degrade the best constraints on ug and Sim, the 
LAI relation is the main source of noise in the measurement of 
cosmological parameters out of the eROSITA cluster sample (see 
Tables vs Table Bl). 

• In turn, combining eROSITA and Planck data gives sensa- 
tional constraints on both the cosmology and the ICM sectors (see 
Table 3). Note that A/^ff^' ~ 7, Acrg = 0.014 and ASl^ = 
0.0039 as well as AaLM = 0.02 and A7LM = 0.07. 

• Other models of primordial non-Gaussianity are analyzed in 
addition to the local shape: the best constraints (with no priors) 
read A/nl — 9, 36, 144 for the local, orthogonal, and equilateral 
shape, respectively, and 7, 19, 115 when Planck priors are applied 
to the other cosmological parameters (see Table 6). 

• Regarding the survey strategy, we find that the all-sky survey 
(/sky ~ 0.658, Tcxp = 1.6 ks) is the optimal choice with respect 
to a possible deeper and smaller survey (/sky ~ 0.15, Tcxp = 7.5 
ks) in terms of parameter estimation (see Table 4). 

• All the results summarized above are obtained considering 
galaxy clusters with mass M > 5 X lO^^h-^MQ (in the fiducial 
model). Lowering this threshold would increase the constraining 
power of the experiments in the Fisher-matrix analysis. This formal 
improvement, however, would be achievable only if robust knowl- 
edge of the X-ray scaling-relations at low masses were available, 
which nowadays is not the case. 

• The forecasted error estimates for the cosmological parame- 
ters are affected by the exact form and redshift evolution of the fidu- 
cial LM scaling relation. Switching from the observationally mo- 
tivated relation by Vikhlinin et al. (2009a) (our reference choice) 
to the one by Mantz et al. (2010a) shrinks the forecasted errors by 
a factor ranging between 10 to 50 per cent, except for /nl^' and 
Sim (whose constraints are barely affected, see Section 11.2). This 
shows that our main results might possibly be conservative. 

• It is impossible to constrain the temperature-mass (TM) scal- 
ing relation on top of our standard 10 parameters (Section 1 1.3). 

• The Limber approximation for the calculation of the angu- 
lar power spectrum is less accurate in the presence of a scale- 
dependent halo bias induced by primordial non-Gaussianity than 
previously established using Gaussian initial conditions. The er- 
roneous use of the Limber approximation on large angular scales 
would optimistically underestimate the marginal constraints on 
/i,°r' by a factor of 4. 

In conclusion, together with the Planck satellite and upcom- 
ing photometric galaxy surveys, eROSITA will substantially con- 
tribute to the simultaneous determination of the cosmological pa- 
rameters and of the X-ray luminosity-mass relation. In particular, 
it will shed new light on the physics of the primordial universe by 
constraining the non-linearity parameter with remarkable accuracy, 
and possibly rule out entire classes of inflationary models. Further 
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studies will focus on the design of dedicated follow-up campaigns 
to further exploit the eROSITA cluster sample. 
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APPENDIX A: ERROR ELLIPSES 

We show here the joint (marginal) 1-a error ellipses for all the pa- 
rameter pairs of our main analysis. Results from different experi- 
ments are color-coded as in Fig. 11, but differently from there re- 
sults are given adopting the external priors on h and Qb, as detailed 
in Section 8.4. Note that the contours obtained from the number- 
count experiment with no redshift information exceed the area in 
the plot in all cases but for the pair fib and h, where they are indis- 
tinguishable from the results of the number counts with photomet- 
ric redshifts. 



APPENDIX B: ADDITIONAL RESULTS 

We report here some additional results which may still be of inter- 
est for the community and shall facilitate comparisons with ana- 
log calculations in the literature. First of all, we show in Table 
Bl the eROSITA constraints on the cosmological parameters that 
would be obtained if the 4 parameters in the LM relation were per- 
fectly known before performing the observations (6-parameter fit). 
These are the most optimistic constraints that can be achieved us- 
ing eROSITA data. Comparing these results with Table 3 shows a 
manifest degradation of the results due to the inclusion of the 4 ad- 
ditional LM parameters: marginal errors on as and f2m are tighter 
by a factor of 10 when the ICM parameters are kept fixed. This 
suggests that extracting cluster subsamples for which precise deter- 
minations of the observable-mass relation could be obtained will 
be a convenient strategy. 

In order to assess the impact of adding /nl'^' to the set of free 
parameters, in Table B2 we repeat the analysis by assuming that 
both the LM relation and /^l*' are perfectly known in advance 
(5-parameter fit). This shows that accounting for primordial non- 
Gaussianity of the local type does not sensibly degrade the con- 
straints on the ACDM parameters, although details depend some- 
what on the individual experiments which are considered. Table 
B2 shows a comparison between eROSITA and Planck, if LM 



and fi^l parameters were perfectly known and where cosmolog- 
ical results for Planck come from a 5-parameter fit of the power 
spectrum of the temperature anisotropics, for Gaussian initial con- 
ditions. Finally, in Table B3 we show how increasing the number 
of free parameters in the fit to the data degrades the marginal con- 
straints on as from the eROSITA cluster number coimts. In this 
exercise, the parameters which are not fitted are assumed to be per- 
fectly known. No prior knowledge of any sort is instead considered 
for the fitted parameters. It is worthwhile mentioning that adding 
the foiu- parameters of the TM relation to our "cosmology-only" 
6-parameter fit gives a marginal error on as which is three times 
smaller than what is obtained by adding the 4 parameters regulat- 
ing the LM relation and its redshift evolution. This occurs also for 
all the other cosmological parameters, although at lower extents. 
Once again this shows that some knowledge of the LM relation is 
the most important piece of information on the ICM physics which 
is necessary to accurately determine the cosmological parameters. 
This holds true also for the measurements of the angular power 
spectrum. 
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Figure Al. Joint I-ct error ellipses for the cosmology sector Color coding is as in Fig. 11. The three top-right panels are a zoomed-in version of the 
corresponding panels in the matrix. 
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Table Bl. As in Table 3 with no priors but marginalizing only over the cosmological sector, fixing the ICM parameters. 
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Figure A2. As in Fig. Al but for the parameters of the LAI scaling relation. The top-right panel is a zoomed-in version of the corresponding panel in the 
matrix. 
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Table B2. As in Table 3 with no priors unless explicitly stated but marginalizing only over 5 parameters of the cosmological sector, fixing both the ICM 
parameters and /J^^^' = 
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Figure A3. As in Fig. Al but for pairs of parameters belonging to different sectors. 
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Table B3. The marginal error on erg as determined from the number-count experiment with photometric redshifts is shown as a function of the number of free 
parameters used to fit the eROSITA data. The last column hsts the parameter showing the largest linear coiTelation coefficient with erg. Note that LM and 
TM stand respectively for the luminosity-mass relation and the temperature-mass relation. The shaded row indicates the reference set-up adopted throughout 
this work. 
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